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In this Issue 

fn 1965. V\W\ Cooley and J\W, Tukey published their landmark "algorithm 
for the machine calculation of complex Fouher series," which qufckty became 
known as the fast Fourier transform, or FFT. The FFT reduced from hours 
to seconds the time required for a digital computer to calculate the frequency 
spectrum of a signal Hewlett-Packard s first FFT analyzer, the HP 5450 A 
Fourier Analyzer, made its debut in 1970, This general-purpose instrument 
could sample input signals between dc and 25 kHz, store the sample values 
in its memory, and later compute frequency spectra or perform other math- 
ematical manipulations of the data 
Over the years, new twists and refinements have made the FFT ever faster, and advancing 
technology has increased processing speeds by a combination of faster digital circuits and parallel 
processing techniques such as pipelining. In this issue you'll find the design stony of a new FFT 
analyzer that uses the latest technology, plus a few of its designers' own ideas, to sample signals 
between dc and 1 00 kHz and compute frequency spectra almost in real time— as many as ten 
spectra per second, fast enough that the user can watch the displayed spectrum change as 
changes occur in the system being analyzed. The HP 3561 A Dynamic Signal Analyzer is designed 
especially for electronic applications such as spectral analysis, network analysis, and waveform 
recording, for mechanical vibration applications such as predictive maintenance, dynamic balanc- 
ing, and spectral analysis, and for acoustic applications such as noise reduction and speech 
signal processing. Our cover illustrates one use of the HP 3561 As near-real-time processing. 
The plot shown is variously known as a spectral map, a stack plot, or a waterfall diagram. It's a 
series of successive spectra displayed simultaneously to show the results of a change of some 
parameter such as the speed of a motor, In this case the system is a fluid-film bearing (like the 
bearings on your cars crankshaft) and the spectral map shows the phenomenon called oil whirl 
becoming oil-whip Instability as the shaft speed increases. For predictive machinery maintenance, 
the HP 3561 A is used just as a physician uses an electrocardiogram. A machines vibration 
spectrum is compared with a normal one taken earlier to reveal changes that might signal an 
impending breakdown because of wear or internal damage. 

The HP 3561 As capabilities are discussed in the article on page 4. On page 6 you'll find 
descriptions of some of its appfications. Its hardware and software designs are described in the 
articles on pages 12 and 17, and the special aspects of its FFT implementation are outlined on 
page 20. The article on page 28 tells how the HP 3561 A's digital filters process so much data 
so quickly. 

Here's an invitation for HP end users and OEMs. Most of our articles are about the design of 
HP products and the contributions made by these products and designs. Wed like to publish an 
occasional article describing how HP hardware and or software contributions have helped you 
make a contribution in your field. We're looking for original technical material that's interesting to 
more than just a few readers. If you have an idea for an article you'd like to write, send an outline 
to Editor, Hewlett-Packard Journal, 3000 Hanover Street, Palo Alto, California 94304 U.S.A. 

-RP. Doian 



What's Ahead 

The design of a new line of general-purpose fiber optic test instruments will be featured in the 
January 1985 issue. The HP 8150A Optical Signal Source produces modulated, calibrated fight 
signals for stimulating fiber optic devices. The HP 8151 A Optical Pulse Power Meter measures 
various parameters of light signais, including peak levels, mesial power level, and extinction ratio. 
A third instrument, the HP 81 51 9A Optical Receiver, is a stand-alone optical-to-electhcal transducer 
that gets users into fiber optic testing at the lowest possible cost. 
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Versatile Instrument Simplifies Dynamic 
Signal Analysis at Low Frequencies 

Analysis of low-frequency signals has many uses in 
electronic design, vibration studies, and acoustic 
measurements. This easy-to-use analyzer covers the range 
from 125 fiHz to 100 kHz and displays the data in several 
useful formats. 



by James S. Epstein 

REAL-TIME ANALYSIS of low-frequency electrical, 
mechanical, or acoustic phenomena can be of great 
value to design and maintenance personnel Until 
the advent of digital signal processing, such analysis was 
difficult if not impossible to dn. As the complexity and 
performance of digital ICs increased, it became possible to 
study low-frequency phenomena in more detail. 

Using fast Fourier transform [FFT T see box on page 20) 
digital analysis techniques, HP introduced several low-fre- 
quency signal analyzers In the 1970s — the HP 5450A 
Fourier Analyzer, the HP 542 OA Digital Signal Analyzer, 
and the HP 3 582 A Spectrum Analyzer Combining many 
of the features of these earlier instruments with new capa- 
bilities made possible by VLSI circuit technology, the HP 
3561 A Dynamic Signal Analyzer (Fig. 1) is the latest 
member of HP's family of FKT analyzers. 

The HP 3561 A is designed to provide high-performance 



dynamic measurement and analysis capability in a small 
portable package. With the increased computational capa- 
bility that exists in today's microprocessors, many features 
can be now incorporated into an instrument that previously 
required a companion controller or computer, Analysis of 
dynamic signals has only recently been available because 
it awaited the development of suitable high-speed digital 
signal processing capability. This capability is incorpo- 
rated in the HP 3 5 61 A in the form of multiple micropro- 
cessors and custom digital integrated circuits. 

Conventional swept -frequency analyzers cannot analyze 
a dynamic signal very well because the signal may be 
changing during the sweep time. In dynamic signal analyz- 
ers, the signal is first digitized quickly and then analyzed 
after the data is captured. The analysis of the data may 
consist of examining the time waveform directly or examin- 
ing the frequency spectrum at evenly spaced intervals or 




Fig. 1. The HP 3561 A Dynamic 
Signal Analyzer combines high- 
performance hardware and ver- 
satile menu-driven software >n a 
portable package for measure- 
ment and analysis of low-fre- 
quency waveforms in electronic, 
acoustic , and vibration applica- 
tions. Besides spectral analysis, 
the HP 3561 A can perform net- 
work measurements and time do- 
main waveform recording. 



4 HEWLETT-PACKARD JOURNAL DECEMBER ^9B4 



)Copr. 1949-1998 Hewlett-Packard Co. 



in Va- or full-octave intervals. In addition, mathematical 
operations can be performed on the data to look at differ- 
ences, ratios, and other computed results. After the data 
has been captured and the desired results computed, the 
next step is presenting the data for display- The HP 3 561 A 
has several display formats, some of which allow up to 
sixty spectra to he viewed simultaneously. 

Key Features 

The HP 3 561 A can measure signals over the frequency 
range from 125 fiHz to 100 kHz. This covers the range of 
slow mechanical and geophysical vibrations to well above 
the 20-kHz audio band. One of the important specifications 
of a dynamic signal analyzer is its ability to resolve signals 
that are closely spaced in frequency. This allows a user to 
detect close-in sidebands on audio signals or vibrations 
caused by two different sources that are close in frequency, 
The HP 3561A can analyze signals with a resolution of B4€ 
^uHz over its entire 100- kHz range. 

Another important consideration in dynamic signal 
analysis is measuring small signals in the presence of large 
signals, such as the distortion of a sine Input. A measure 
of an analyzer's ability to do this is called dynamic range; 
it specifies the maximum difference in signal levels that 
can be measured and guarantees that any internally gener- 
ated distortion and spurious products will be below the 
minimum measurable level. In the HP 3561 A the dynamic 
range is 80 dB t which means that a primary signal and a 
secondary signal une ten-thousandth of the amplitude of 
the primary signal t:an be measured simultaneously. The 
accuracy of the measurement is also important. Using an 
internal calibration process, the HP 3561 A has an 
amplitude accuracy of ±0.15 dB over the entire HP class 
B environmental range (0 to 55*C, <95% relative humidity 
at 40°C, and altitudes <4570 m). 

Applications 

The applications for a dynamic signal analyzer are very 
diverse- However, the HP 3561A Dynamic Signal Analyzer 
is specially designed to cover three areas of application. 

The first area is general electronic measure men Is in the 
audio frequency range. One common measurement is I lie 
distortion of a pure tone. The HP 3 561 A can measure dis- 
hirlinn prnihn Is I hat are 80 dB below Ihe fundamental and 
has a built-in marker function thai computes the total j i 
tort ion of up to 20 harmonics when the marker is placed 
on the fundamental. Another electronic measurement is 
treqmmcy response. The HP 35 til A has an internal source 
that, combined with the internal trace arithmetic, can he 
used to measure both amplitude and phase as functions of 
frequency. General spectral analysis is another measure- 
ment dime in the audio band to determine the exact fre- 
quency and amplitude of an input signal or to determine 
ihe frequency and amplitude of a spurious input. Transient 
analysis is another measurement often needed in the elec- 
tronic and mechanical areas. The 40, 960- word buffer in 
the HP 3561 A allows 0.8 second worth of data to be cap- 
lured for signals in the (Mo-20-kllz frequency range. 

Another major application area for dynamic signal 
analysis is measuring mechanical vibration using a trans- 
ducer that converts either displacement, velocity, or accel- 



eration into an electrical signal that can be measured by 
the analyzer. The analysis of vibration can be extremely 
useful in the maintenance of machinery. By examining the 
vibration in the frequency domain, small vibrations caused 
by worn-out parts can be detected even in the presence of 
large vibrations that are normal to the operation. In this 
manner, the machine can be shut down for repair only 
when necessary, saving costly down time, hi addition to 
detecting when a problem has occurred, potential causes 
can be identified. (See box on page 6 for an example), 

A third application area for which the HP 3561A is par- 
ticularly well suited is acoustic measurements. Measuring 
acoustic signals requires special techniques. The human 
ear has its owm frequency response, which must be con- 
sidered when making measurements. An internal A- 
weighting filter allows the analyzer to duplicate the weight- 
ing applied by the human ear. thereby allowing the noise 
levels measured to be compared to their levels as perceived 
by a human observer. Another important consideration in 
acoustic analysis is frequency resolution. The human ear 
i annot detect iwo frequencies separated by a fixed differ- 
ence at high frequencies as well as it can at low frequencies* 
Its response is logarithmic and is best characterized by 
using ^i-octave analysis bands that become wider as the 
center frequency increases. These analysis bands have been 
used for a long time to measure acoustic signals and many 
noise specifications are written in terms of Va-octave inter- 
vals as defined by ANSI specifications. The HP 3561A can 
display its results in Va-octave. full-octave, or linear fre- 
quency increments to allow the user to determine which 
mode is most useful for a particular application. 

Friendliness 

With all of the capability that can be packed into an 
instrument such as the HP 3561 A, It could be difficult to 
use these features if the instrument were not specifically 
designed for ease of use. A measurement should be easy 
to make without requiring a lot of setup to get accurate and 
quick results. This is made possible in the HP 3561 A by 
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Dynamic Signal Analysis 
for Machinery Maintenance 



Using the HP 3561 A Dynamic Signal Analyzer to analyze the 
vibrational signature of machinery can allow identification of prob- 
lem spots before ihey cause machinery shutdown. In the time 
domain, vibrations from many sources in a machine can add to 
create compficated waveforms. When viewed m the frequency 
domain, however, individual problems such as imbalances, gear- 
mesh anomalies, or bearing defects can usually be identified 
(Fig. 1). By saving the spectra of a healthy machine in the HP 



a damaged gear as shown by the excessive response at the 

gear-mesh frequencies (Fig. 2b). 

In addition, the HP 3561 A has a spectral map display format 
that is particularly useful in servicing rotating machinery These 
spectral maps can be created by adding a third dimension to 
spectral measurements, such as time, r/min, or load The HP 
3561 A's map mode can be used to identify changes or trends 
in the spectra of machines, circuits, or other devices. In the exam- 
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Fig, 2. faj Baseline spectrum for a machine, (b) Later spectrum for machine of (a) showing 
the excessive response at the gear-mesh frequencies caused by a damaged gear. 



3561 A r s optional one-megabit bubble memory, a set of baseline 
measurements is formed that can be used for later comparison. 

Fig 2a shows an example of a normal baseline measurement 
for a machine The same machine at a later time has developed 



pie shown in Fig 5 on page 10. a collection of measurements 
made during a machine runup forms a spectral map, illustrating 
changes in machine vibration with an increase in rotational 
speed. 
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several features coupled with softkey menus. An aut orange 
algorithm selects the correct range to place the input signal 
within 4 dB of full scale, Autocalibration is periodically 
performed so that the measurement accuracy is maintained 
within the _D.15-dB specification. The first measurement 
covers the entire 100-kHz frequency range so thai all fre- 
quencies can be examined. Then the user can zoom in on 
a particular range of frequencies by using the HP 3561 A' s 
marker function to identify the center frequency and press- 
ing the DEFINE CENTER softkey in the menu obtained by 
pressing the FREGuency key on the front panel. After the 
center frequency is selected, the desired analysis band can 
be selected in the FREQ softkey menu by pressing the DEFINE 
SPAN soil I 

Several special marker functions are provided to make 
the analysis of the measurement easier. The relative marker 
function allows any signal to be selected as the reference 
for either amplitude or frequency or both. When this marker 
is used, all measurements are scaled relative to the selected 
reference. This function is particularly useful when 
measuring how far the amplitude of a spurious signal is 
below that of the desired signal. 

Another special marker function will, mark up to twenty 
harmonics on the display screen with a special grid after 
the fundamental frequency is entered by the user [Fig, 2). 
In addition to marking the harmonic components on the 
screen t the total harmonic distortion is computed for the 
harmonics marked, A similar marker function is provided 
tor marking Sidebands and computing the relative power 
in the sidebands with respect to the carrier, 

A very common type of analysis is computing the power 
in a specified frequency band. This is done w T hen analyzing 
noise power in different frequency regions. The band- 
power marker function allows the user to position right 
and left markers to denote the band of interest and have 



the power in the band computed. 

Saving Instrument States 

Some specialized measurements may require many of 
the HP 3561 A*s special features to be selected. When this 
is the case, the storage and recall of instrument states in 
nonvolatile memory can he a very handy feature. Once the 
instrument is set up as desired, its state can be stored in 
one of six locations. If the optional bubble memory' is in- 
stalled, up to 127 instrument states can be saved, At a later 
time, say at a remote measurement site, the desired state 
can be quickly recalled when it is time to make the measure- 
ment. This is not only convenient, but can help prevent 
errors in setting up the instrument. 

Narrowband Measurements 

The narrowband measurement is the basic analysis mode 
nf the HP 3561 A. This mode provides a frequency spectrum 
of 400 spectral lines equally spaced over the chosen 
analysis band. The analysis band can be centered anywhere 
in the instrument s frequency range from dc to 100 kHz 
with 0,25-Hz resolution and the frequency span for a 
zoomed measurement can be as small as 256 mHz. This 
gives a resolution of 640 pHz anywhere in the spectrum 
up to 100 kHz. In addition, when making a baseband mea- 
surement [not zooming around a center frequency), the 
analysis band, or frequency span, can be reduced to 10.24 
mHz. which gives a maximum resolution of 25.6 ixHz. 
Speed. Speed was an important consideration in designing 
the hardware to make narrowband measurements. When 
making adjustments to electronic equipment, such as re- 
ducing distortion components in a device under lest f ihe 
display must be updated fast enough to give useful feedback 
to the person doing the adjustments, When averaging many 
spectra, the total measurement time is greatly dependent 
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upon the time required to generate each spectrum. The HP 
3 561 A can update a spectral display 10 times per second 
and can average up to 20 spectra per second. See the article 
on page 12 for design details about how this performance 
was achieved. 

Windows. The shape of the analysis band is determined 
by the window applied to the data. In dynamic signal 
analysis, the band shape is selected by windowing — mul- 
tiplying the data in the time domain by a weighting func- 
tion, This is equivalent to filtering the data in the frequency 
domain. Just as it is useful to have different passband 
shapes in a conventional sweptTrequency analyzer, it is 
also useful to have different windows in a dynamic signal 
analyzer. The HP 3 561 A offers a choice of four window 
functions — flattop, uniform, Manning, and exponenlial. 
The fiattop window gives maximum amplitude accuracy 
at the expense of frequency resolution. Its bandwidlb is 
approximately four times that of the uniform window, 
which weights all of the time samples equally. The uniform 



window is used for signals that are periodic or are fully 
contained within the analysis band (such as a Iransient or 
an output of one of the HP 3561A*s internal noise sources] . 
and it provides maximum frequency resolution. 

The Hanning window provides a good compromise be- 
tween frequency resolution and amplitude accuracy. The 
window flatness is within 1,5 dB, compared to 0,01 dB for 
the flattop window. Its bandwidth, however, is only 1-5 
limes the maximum available resolution. 

The exponential window is useful when measuring 
transients. Coupled with the HP 3561 A 's internal impulse 
noise source, this window allows the HP 3561A to make 
impulse response measurements. The shape of the expo- 
nential window can be adjusted to fit the data that has 
been collected. 

Measuring Phase. Another distinction between dynamic 
signal analyzers and conventional swept-frequencv analyz- 
ers is the ability of the former to measure phase. There are 
also applications for measuring phase with respect to a 
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trigger, such as when balancing rotating machinery. Here, 
the measured phase is used to determine where In the 
rotation cycle an out- of- balance condition occurs. By using 
the HP 3 561 As internal calibrator, the measured phase 
value is corrected for an accuracy of ±2 degrees up to 10 
kHz and ±10 degrees from 10 kHz to 100 kHz. 
Averaging, Averaging can be used to increase the useful- 
ness of £he narrowband measurement. Several types of av- 
eraging are available. When measuring the spectral charac- 
teristics of random data, the accuracy of the results is sub- 
ject to the statistical properties of the data. One method of 
reducing the variance of the results is to perform rms or 
power averaging of the data. Another type of averaging that 
is useful for reducing the amount of noise on the signal is 
time averaging. This type of averaging requires a synchroni- 
zation signal, which can be either an external trigger or a 
part of the signal being measured. Two other types of av- 
eraging provided by the HP 3561 A are rms exponential and 
peak hold. The first is a type of rms averaging that weights 
earlier averages less than the current measurement. This 
allows a continuous update of the display while averaging. 
The second is used when the maximum value in each fre- 
quency bin is desired for a series of measurements. 
Triggering, Triggering is important when measuring phase 
and in time averaging. Therefore, the HP 3 561 A has several 
triggering modes. The source of the trigger can be the input 
signal, an external TTL signal, the internal source, or the 
internal trigger (free run). In addition, the instrument can 
be armed for a single trigger or armed to respond each time 
a trigger is received. Pretriggermg of up to 8 records in 
narrowband mode or 40 records in time capture mode (dis- 
cussed later) is provided. Post- triggering of up to 1023 rec- 
ords can be selected, Pre- and post-triggering are very useful 
when a time difference exists between a known event ;uid 
the data that is desired. When testing loudspeakers, for 
example, a time delay exists between the time the impulse 
excitation is generated by the instrument and when the 
impulse response is detected at the speaker. 
Simultaneous Display. The HP 3561A has the ability to 
display a full 100-kHz baseband spectrum and a zoom spec- 
trum simultaneously [Fig. 3), This capability is provided 
by the use of custom digital filters, which are discussed in 
the article on page 28. One way that this display format 
can be used is to have the full 100-kHz spectrum on the 
upper half of the display and the zoom spectrum on the 
bottom. Moving the marker to the desired signal on the 
upper display and pressing the CTR FREQ softkey in the 
marker menu will move the center frequency of the zoom 
spectrum to that of the desired signal. This can be per- 
formed repeatedly if desired to zoom in on several different 
signals for closer examination. 

Time Capture Measurements 

Time waveform capture can be used for transient analysis 
and the analysis of nonstationary signals such as speech. 
It also can be used to increase the real-time rate of the 
measurement to 100 kHz until the buffer is filled. 

The major difference between narrowband and time cap- 
ture measurements is that in time capture mode the entire 
lime record is captured and stored in memory before any 
processing is done, If the bubble memory option is in 



stalled, the time record can be retained even while the 
instrument is turned off. The time record can be defined 
to be from 1 to 40 records of 1024 samples each. The full 
triggering capability of the HP 3 561 A can be used with 
time capture mode so that up to 40 records of pretrigger 
data can be taken, allowing the examination of data before 
the occurrence of a particular event. The ability to preserve 
the original time data before processing allows the center 
frequency and span for zoom analysis to be selected after 
collecting the data. This can not only save a lot of data 
collection time, but may be the only way to collect data 
when the input signal is of a transient nature. 

In time capture mode, the data is first collected within 
the frequency range specified for collection. If the data is 
collected in baseband mode T it can be translated to any 
center frequency up to the maximum frequency collected 
and zoomed up to a factor equal to the number of records 
collected. 

The postcolh M Mil translation can occur as many times 
as desired. In addition, the type of window applied to the 
data can also be changed to allow trade-offs between fre- 
quency resolution and accuracy. The postcol lection pro- 
cessing is done by selecting the time capture data memory 
as the input to the digital filters. Because the data has 
already been collected, it can be run through the digital 
filters at maximum speed. Different start points in the data 
can be selected so that only desired data is analyzed and 
changes in spectral parameters can be examined at different 
points in the record. In addition, an average magnitude 
display can be selected to examine the average of all I he 
data in the time capture memory at one time. 

When the time capture mode is combined with other HP 
3561A features, very useful measurements can be made. 
For example, the average magnitude display can be used 
in I he upper half of the display to show the full frequency 
range of the captured data, and the lower half of the display 
can then be set up to scroll through the captured data with 
or without zoom. When a user desires to look at the entire 
time buffer at one time, the spectral map display mode 
(discussed later) can be used. Whether scrolling through 
the data or mapping it, an increment can be set up to specify 
the percent of the dat<i m "nl lo be incremented for the 
next display. In this manner, a spectral map can be dis- 
played so that each new trace in the map is offset by one 
to 1024 data points. This can be important when examining 
the decay of selected frequencies. 

Third- and Full-Octave Analysis 

For certain applications* most notably in the acoustic 
field, proportional frequency resolution is desirable. In this 
type uf spectrum analysis, the signal is analyzed by a bank 
of fillers with logarithmically spaced center frequencies 
and band widths proportional to the center frequency, The 
most widely used analysis of this type is at '/Voctave and 
full -octave intervals, The HP 3561 A computes ihe Woctave 
jEid full-octave spectra digitally by combining narrowband 
frequency data tfiken on ihree separate analysis spans. To 
do this, the custom digital filters are set up to collect data 
at ihree sample rates: f tt , f s /10. and f 9 /100. The three data 
buffers are then transformed using the FFT. and a weighted 
sum of adjacent FFT filter bins is computed to form each 
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1 <-octave or iuJl-«ji lave Biter bin. 

The frequency response of the octave filter shapes is 
determined by the windowing function used, Ihe number 
of FFT bins summed, and the weighting function applied 
to the FFT bins. In the HP 3561 A, the Hanning window 
function is used for all of the Va-octave and full-octave 
filter shapes. The number of bins summed varies with the 
center frequency of the octave filter. The weights applied 
to the adjacent FFT bins equal unity except for the first 
and last bins. 

The octave filter shapes are designed to meet ANSI and 
IEC standards, Because of the number of different octave 
spans allowed in the HP 35filA t 168 special filter shapes 
had to be designed and verified. Since this would have 
been a formidable task if it had been done manually, a 
computer program was written to find the filter weights. 
This program optimized the filter weights based on the 
noise equivalent bandwidth and the geometric mean center 
frequency specifications, which were the severest design 
constraints, A closed-form solution was used for the noise 
equivalent bandwidth and the geometric mean center fre- 
quency was optimized by using a binary search. The pro- 
gram completed the design task altera month of unattended 
operation on an HP 9000 Series 200 Computer! 



available options. The major functional areas are Measure- 
ment, Trigger, HP-IB, Display, Marker, Instrument State, 
and Input. These groups are clearly indicated on the front 
panel so that the user knows which keys to press for any 
desired function. In addition to physical grouping the dif- 
ferent types of keys are color coded so that related turn I i n 
are closely associated. Within each functional group are 
specific keys that either perform a specific action or bring 
up softkey menus for further selections (Fig, 4). The softkey 
menus displayed along the right edge of the CRT screen 
show only the choices that are available for the mode 
selected. 

The keyboard is made up of a new conductive rubber 
key set that is custom molded for the HP 3561A's front 
panel This keyboard allows multiple colors and very sharp 
lettering to be molded into the keys. The keys have a tactile 
feedback that Indicates that a key has been pressed. In 
addition to their attractive appearance, the conductive rub- 
ber keys provide a ruggedness that is commensurate with 
the portability of the instrument These rubber keys are 
very difficult to damage, even If objects are dropped on 
them. Furthermore, they provide resistance to dust and 
moisture and are very easy to install, reducing the assemhl \ 
time and the cost of the instrument, 



Front Panel 

The front panel on the HP 35 til A is specifically designed 
for ease of use and ruggedness. The user interface with the 
instrument is through keys grouped by function and as- 
sociated softkey menus that direct the user through the 



Softkeys 

The softkey menu interface allows the user to be guided 
through the selections required for a setup. When a key, 
such as the MODE key in the Measurement group, is pressed, 
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Fig. 5. Spectral map display. Up to 60 magnitude traces can be displayed simultaneously for 

studies of changes tn the spectrum as a function of ttme, r/mtn r etc. These maps can be saved 
and recalled tater from the optional nonvolatile bubble memory, 
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the associated choices then available to the user are dis- 
played along the right edge of the display (see Fig. 1). For 
the MODE key, the resulting choices are NARROW BAND, 
THIRD OCTAVE, FULL OCTAVE, TIME CAPTURE, EXT SAMP ON 
OFF, DERNE PULS/REV, and TEST SELECT, This is a good 
representation of the different types of softkeys available. 
The first four softkeys make up a group of selections, only 
one of which is allowed at a time. When the selection is 
made, the associated softkey label is intensified so that it 
is obvious which softkey is presently in effect. The next 
type of softkey is the on/off softkey or two-state softkey 
This softkey allows the user the choice of having the labeled 
function either on or off. For example, the EXT SAMP ON OFF 
softkey allows the user to supply an external sample clock 
instead of using the internal one. When an on off softkey 
is pressed, it toggles between on and off with its current 
state indicated by either the ON or OFF in its label being 
intensified. Another type of softkey is typified by the DE- 
FINE PULS REV softkey. This softkey requires a numeric 
entry, as do all softkeys that have the word DEFINE in their 
labels. When the label for this type of softkey is intensified, 
any numeric entries made will be for the parameter indi- 
cated on its labeL If more than one DEFINE softkey label is 
available on a menu, only one will be intensified at a time. 
The units for the entry will be displayed on the softkey 
menu af lei any numeric key p| the front panel's keypad or 
the intensified DEFINE softkey is pressed. 

The final type of softkey available is the select softkey 
as represented by the TEST SELECT softkey. When a softkey 
with SELECT in its label is pressed, another menu is dis- 
played with more selections. (See article on page 17 for 
more detail about the HP 3561 A's software.) 

Display Formats 

Once the data is collected Y the next step is to display it 
in a manner that best allows the user to analyze the results. 
Bet :a use of the many different types of measurements pos- 
sible, ihe type c j I analysis performed will vary depending 
upon the application. Fur this reason , several different dis- 
play modes are provided. We have already discussed the 
upper/lower display formal in which two separate traces 
of any type can be displayed simultaneously, one above 
the other. There are many uses for this and il is the format 
that the instrument powers up in. Some examples are look- 
ing at magnitude and phase at the same time, observing 
the time waveform along with its magnitude or phase spec- 
trum, or viewing the full 100- kHz basehand spectrum along 
with a zoom spectrum, In the time capture mode, viewing 
a compressed time record and a magnitude spectrum to- 
gether allows the user to see which part of the large time 
record is current I \ being processed. 

A single display format is also provided so that a display 
can be expanded to lake up the entire screen when only 
one trace is necessary to display the measurement. Another 
display format that is very useful is the front/hack display. 
This display also allows any two traces to be viewed simul- 
taneously, However, they are both expanded to fill the 
entire screen with one positioned behind the other. The 
back trace is 1m 1 1 the brightness of the front trace so that 
the two traces can be easily compared. The front and back 
traces can be switched by pressing the NEXT TRACE key 



The most versatile display format is the spectral map 
display (Fig, 5) obtained by pressing the MAG MAP CLEAR 
softkey in the display format menu. Up to 60 magnitude 
traces can be simultaneously displayed, one behind the 
other, with a hidden-line algorithm applied so that the dis- 
play is more viewable. The traces can be positioned either 
directly behind each other or slightly offset along the hori- 
zontal axis so that data that may be blocked by the hidden- 
line algorithm can be seen, Through the use of the NEXT 
TRACE key and the Marker group keys, the full accuracy of 
the instrument is available while viewing any one of the 
60 traces. The NEXT TRACE key allows stepping through the 
map with the selected trace intensified. The marker func- 
tion then reads out the x and y values for the location of 
the marker to full precision. If it is desirable to expand the 
selected trace to fill the screen, a toggle key allows the user 
to switch between the spectral map and the selected trace. 
The data in the map is preserved in the HP 35GlA's memory 
until a new map is generated by pressing the MAG MAP 
CLEAR softkey again- If the bubble memory option is in- 
stalled, the entire map can be stored in nonvolatile m^ 
for recall later. 

Portability 

Dynamic signal ana lysis applications are not always con- 
fined to the iab bench. These measurements are often made 
in remote locations such as on aircraft or at remote com- 
munications sites. Vibration analysis can be performed on 
machines that are located on shop floors, in power plants, 
in chemical factories, or in submarines. Hence, the HP 
3561 A is packaged much the same as a portable oscillo- 
scope, it weighs 33 pounds. A convenient carrying pouch 
is attached to the top of the instrument to store cables, 
manuals, and other equipment. 

Portability implies more than just a package; If a lot of 
additional equipment is required, portability is decreased. 
For example, when the instrument is carried to a measure- 
ment site, it may be desirable to bring the data back to a 
central computer for further analysis or storage. The op- 
tional bubble memory for the HP 3561 A allows up to 127 
traces to be stored and retained even when the power is 
turned off. The traces can be recalled at any time for further 
analysis or transfer to another computer. 

When making vibration measurements, some kind of 
transducer is necessary. To reduce the number of transduc- 
ers that have to be carried around with the instrument, 
inexpensive integrated circuit piezoelectric fICP) ac- 
celeTometers can be left attached to the machine under t«-s( 
The current source required to power them is built into 
the HP 3561 A as a standard feature. 
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Hardware Design for a Dynamic Signal 
Analyzer 

by James S. Epstein, Glenn R, Engel, Donald R. Hiiler, Glen L. Purely, Jr., Bryan C. Hoog, 
and Eric J. Wicklund 



■^ ESIGNING THE HARDWARE for the HP 3561 A 
I Dynamic Signal Analyzer involved the constraints 
b^r of portability, adequate shielding of the signal pro- 
cessing circuits against electrical noise, a compact power 
supply, high-speed digital circuitry* and an accurate clear 
display of data and calculated results, Supporting and 
simplifying some portions of this hardware design are an 
extensive softkey menu structure and a sophisticated 
operating system (see article on page 17). 

Hardware Design 

A block diagram of the HP 3561 A is shown in Fig, 1. 
The front end performs the signal conditioning required 
to scale and filter the data so that it can be digitized by the 
analog-to-digital converter (ADC). After the signal is dig- 
itized, it is further processed by the digital filters to provide 
frequency translation and increased resolution. The main 
processor controls the entire instrument to perform the 
tasks required and communicates with the display to show 
the results of the measurement. The power supply provides 
the power to run all the internal circuitry and the non- 
volatile memory stores traces, setups, and time records for 
further analysis. An optional bubble memory can be added 
to extend the available memory. 

Front End 
The input amplifier circuitry [Fig. 2] in the HP 3561 A 



provides a floating, shielded input connection for rejection 
of common mode and ground loop induced signals, A 
grounding switch is included for chassis-referenced mea- 
surements. The signal path does not become ground refer- 
enced until it leaves the ADC and crosses an isolation trans- 
former. Therefore, all circuitry through the ADC runs on 
separate floating power supplies. 

The input signal passes through a 0- P 20-, or 40-dB at- 
tenuator (depending on range) into a JFET input amplifier 
stage. Four fixed-gain amplifiers, each preceded by a pro- 
grammable attenuator, provide a total gain of 38 to -40 dB 
in 1-dB steps. This gives the HP 3561 A a range of full-scale 
input values from 27 to -51 dBVrms in 1 -dBVrms steps. 

The full-scale rms output voltage level of the input 
amplifier is 228 mV. Gain is distributed through the 
amplifier chain to maintain a reasonable signal level lor 
best noise figure and distortion performance. The resulting 
noise floor specification is -141 dBV/V Hz for frequencies 
less than 1 kHz and - 150 dBV/VHz above 1 kHz. Distortion 
is less than —80 dB relative to full scale. Internal offset 
voltages are corrected by an eight-bit DAC-controlled auto- 
zero circuit, which sums an opposing offset voltage into the 
last amplifier stage. A binary search algorithm is used by 
the firmware to set the D-to-A converter, 

Following the input amplifier is a passive nine-pole anti- 
aliasing filter w T ith a rolloff above 100 kHz. The signal is 
then buffered, sampled, and digitized by the ADC. 
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Fig, 1. Hardware block diagram of the HP 3561 A. 
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Additional circuitry includes the 4-mA current source 
for use with external integrated circuit piezoelectric (ICP) 
aecelerometers, An internal A- weigh ting filter can be selected 
to provide an ANSI standard simulated "human ear" transfer 
function in the amplifier path, 

A-to-D Converter 

Requirements of the .ADC are set by the bandwidth and 
dynamic range specifications of the instrument. The anti- 
aliasing filter technology used allows band widths up to 
approximately 0.4 times the sample frequency, which is 
fairly close to the theoretical Nyquist limit of 0.5. This 
results in a sample rate of 256,000 samples per second for 
the instrument bandwidth of 100 kHz. Dynamic range con- 
siderations imply a resolution of 13 bits and maximum 
non linearity of T 4 LSB (least-significant bit). 

The HP 3561 A*s ADC achieves its high accuracy by using 
a technique known as residue conversion. The final result 
is formed by performing successive passes or partial A-to-D 
conversions of increasing accuracy and decreasing ran^e. 
Each new pass starts by converting the ADC output (the 
result of the previous passes) to an analog signal using a 
DAC and then subtracting it from the ADC input. The result- 
ing residue represents the errors in the previous approxima- 
tion. The residue is then converted with a resolution and 
accuracy much greater than for the previous pass, The re- 
sult is added to the previous approximation of the ADC 
input to produce a much more accurate approximation, 

A partial block diagram of the ADC is shown in Fig. t 
This two-pass converter uses the same local ADC fan ADC 
within an ADC) to perform the conversions for each pass. 
During the first pass the ADC input signal is routed to the 
local ADC through an 8:1 attenuator. The result of the local 
ADC conversion is the first -pass approximation of the input 
and is stored in latch LI, This first-pass approximation is 
the input to the upper etghl bits of tht- 13-bit DAC, which 
converts with a linearity equivalent to V* LSB of 13 bits. 

The second pass begins as the first-pass approximation 



is latched and the DAC begins to settle- The DAC output 
is subtracted from the ADC input signal and the result is 
the residue representing the errors in the first-pass approx- 
imation. The residue is multiplied by four by an amplifier 
stage before the local ADC performs the second -pass ap- 
proximation. The gain ratio between the first and second 
pass is 4/(1/8] or 32. so the second pass resolves with 32 
times greater accuracy, but has a range 32 times less, The 
second-pass approximation is then added to the first-pass 
approximation in adder Al, but is shifted down five bits 
to account for the gain ratio of 32 between the first pass 
and the second pass. This final result is stored in latch L2 
before being transmitted to the digital filters. 

The most important aspect of a residue converter is the 
cancellation of first-pass errors. If an error E is made in the 
first pass, the DAC is loaded with the correct approximation 
to the input plus E. This is subtracted from the input so 
that the residue is only -E. This is then converted and 
added to the first- pass approximation to produce the final 
result so that — E, the residue, cancels the E term contrib- 
uted by the first-pass approximation. 

An important parameter in a residue converter is the 
overlap between passes. The HP 3561A 1 s ADC has three 
bits of overlap between the first and second passes, mean- 
ing that the upper three bits of the second-pass approxima- 
tion have the same binary weight as the lower three bits 
of the first-pass approximation. The overlap is a measure 
of the magnitude of residue that can be converted by the 
second pass wiihoul overloading. If the first-pass approx- 
imation were perfect, then only quantization errors would 
be present in the lirsl-pass approximation and no overlap 
would be required. 

The significant sources of linearity errors in this design 
are the DAC errors, the ADC errors on the second pass, and 
the analog gain matching between the first and second pass 
(how exactly the gain ratio of 32 is held). The gain matching 
determines how much the first-pass errors contribute to 
the nonlinearitv of the converter. Since, in the best case. 
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Fig, 2. input circuitry of (he HP 3561 A 
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Fig. 3. Partial block diagram of the ADC used in the HP 
3561 As front end. 

the first-pass approximation still has quantization errors 
[the conversion is only eight bits, not 13), the gain matching 
is always important. 

Most ADCs (and DACs] have particularly nasty forms of 
non linearity. Common forms are step and impulse voltage 
transfer function errors. A step error is a constant error thai 
is only added to the output when the input is m a certain 
range. This can occur when I he analog weight associated 
with a particular bit is slightly in error. A common occur- 
rence is when the MSB [most-significant bit) of a DAC is 
in error. In this case a constant error is present wmenever 
the input is greater than half scale, An impulse error is an 
error that is very localized. These nonlinearities result in 
distortion products that do not decrease [and can even rise) 
as the input level is lowered. This type of behavior is not 
desirable in a spectrum analyzer. One technique particu- 
larly effective for reducing this behavior is dithering (also 
discussed in the article on page 28). 

Dithering is the addition of an extra signal to the input 
to "smear out 11 nonlinearities. This causes any particular 
input level to correspond to different points on the voltage 
transfer function at different times. In the long term the 
net transfer function at a particular input level becomes 
the weighted average over the neighborhood of that level 
in the original transfer function. The weighting function 
is the probability distribution function for the dither signal. 
Another way of looking at the effect of dithering is that 
harmonic products tend to be converted into intermodule 
tiun products between the input and dither signals when 
the nonlinearities are localized as typical ADC and DAC 
nonlinearities are. This intermodulation effect implies that 
the dither signal should be broadband enough to distribute 
these products into many frequency bins so that the abso- 
lute magnitude at any frequency is much lower than if the 
products were distributed over only a small group of fre- 
quencies. This is equivalent to randomizing the errors so 
that they are not correlated with the input signal. 

A natural choice for a broadband dither signal is random 
noise. However, dithering only smooths out nonlinearities 
over a range the size of the dither signal, so distention 
improvement is negligible for input levels greater than the 
dither level. This implies that \ht< dither signal must be of 
a greater amplitude than the input signal to receive any 



distortion improvement. One way to increase the level of 
dither while not increasing the noise floor is to subtract 
the dither signal after the conversion. However, this rules 
out many sources of dither signals since the dither signal 
must be easily available in digital form, Hence, the use of 
pseudorandom noise is a natural choice. 

The HP 3 561 A 's ADC uses this form of dithering in a 
special way. For the dithering signal to he available both 
as an analog signal and as a digital signal implies that some 
sort of additional data conversion must be present, usually 
a DAC, It is possible to exploit the residue converter topol- 
ogy and achieve the same result without requiring another 
DAC. Pig, 4 shows the ADC with the dithering circuitry 
added. 

Since the first- to-second-pass gain match is not perfect, 
it is necessary to dither the first -pass conversion to ran- 
domize first-pass non linearity errors. Here the dithering 
signal Nf looks just like any other first-pass error and so 
is subtracted out in the combination of the first-pass and 
second-pass conversions. Since this dither signal appears 
in the DAC input and hence the residue, the DAC and 
second-pass conversion are somewhat dithered. However, 
since the first pass quantizes Nf, the dither signal has only 
several values and is not able to smear out nonlinearities 
adequately. A more continuous range of dither values is 
provided by adding pseudorandom noise to the first-pass 
conversion using adder A2. This dither signal also acts like 
a first-pass error since it is added before latch Ll . The result 
is that all of the sources of nonlinearity are dithered without 
the use of an extra DAC or ADC. The only extra circuitry 
required is a simple analog noise source Nf (actually one 
bit of pseudorandom noise passed through an RC filter to 
prod nee a more uniform probability distribution), the 
pseudorandom noise source, and adder A2, which per- 
forms only carries since the noise source has a zero mean. 
The circuitry within the shaded area in Fig. 4 is realized 
with a gate array in the HP 3561 A, keeping the added cost 
of this dithering scheme very low. 

Processor 

The main processor of the HP 35 61 A is an 8-MHz 68Q00 
microprocessor with 393,216 bytes of program ROM. This 
processor controls the hardw F axe and the user interface and 
performs the display processing and much of the data pro- 
cessing (see Fig. 1). It interfaces to the other hardware sub- 
systems via two different 16-bit buses: the processor bus 
and the RAM bus. The processor bus interfaces the main 
processor to the standard CMOS and optional bubble nan- 
volatile memories t the display hardware, the local oscil- 
lator [LO) f the HF-JB (IEEE 488) interface, the source, and 
the keyboard, 

The HP 3561 A has 1 3 1 ,072 bytes of dynamic RAM, which 
is accessed through the RAM bus. This RAM is organized 
as B4K 18-bit words, but its data can be accessed as bytes 
by the miiin processor. This RAM is used for the processor's 
stack and program variables and for all the data buffers 
needed during the data processing. The RAM bus is con- 
trolled by a bus arbiter and the dynamic RAM controller. It 
has DMA ports to the dedicated FFT processor and the 
outputs of the digital filter control chip. The digital filter 
control chip performs DMA accesses using two address 
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counters which control where in RAM the data is written. 
The RAM bus also has some control and status ports used 
by the main processor for interfacing to the FFT processor. 
the digital filler control chip, the digital filter DMA address 
counters, timing logic, and the interface to the HP 3561 A *$ 
front-end circuitry and ADC, 

Nonvolatile Memory 

The standard nonvolatile memory in the HP 3 561 A con- 
sists of 8K bytes of CMOS RAM. This RAM is used to store 
data for two traces and six instrument setups, with one of 
the setups being the last state before power-down. The 
CMOS RAM is backed up by a three- volt lithium battery, 
which allows it to retain data for up to five years. 

There is also a bubble memory option that can add a 
megabit of nonvolatile memory. This allows users to store 
extended time records, as well as more setups and traces, 
The bubble memory can be accessed by a simple processor 
interface and can transfer data at a rate of 12.5K bytes/s. 
These speeds allow the user to retrieve an extended time 
trace of forty records in 0.7 second, worst case. 

Noise Source 

The noise source in the HP 3 561 A can generate three 
types of signals: periodic, impulse, and random. All three 
of the signals are band-limited and band-translated. Like 
the tracking oscillator in a swept analyzer, Ihe noise source 
always tracks the selected frequency span. 

The impulse source has the added capability of being 
triggered from the HP 3561 A's front panel or via the HP-IB. 
and it can skip time records. These features are useful when 
transient inputs are required to settle before a new measu re- 
nunt is made. The impulse source can skip up to 32,766 
records before generating another impulse. 

Digital Filter 

The digital filter section is based on custom ICs de- 
veloped especially by HP for dynamic signal analysis. Con- 
trolled by the main processor, this section tillers the dig- 
itized input waveform to form the appropriate data records 



for the FFT processor. See the article on page 28 for a 
detailed discussion of the design constraints and realiza- 
tion of this section. 

Display 

The HP 3561 A uses a raster-scanned magnetic-deflection 
CRT for display of measurement results* This Is a departure 
from some other Hewlett-Packard signal analyzers, which 
use a vector-scanned electrostatic-deflection CRT display. 
The raster-scan technology offers significant advantages for 
the HP 3561 A. Probably the most important is compact 
size, With magnetic-deflection CRTs, large electron beam 
deflection angles are easily achieved. Because of this, mag- 
netic-deflection CRTs can be much shorter than eleclrasta- 
t j c-def lection CRTs for a given viewing area. The space 
saved by using a magnetic-deflection display helped fit the 
HP 3561 A 's extensive measurement capability into an os- 
cilloscope-sized package. Another advantage of the raster- 
scan technology is low cost. The HP 3561 A display, includ- 
ing drive electronics, costs about as much as an electrostatic 
CRT by itself. 

Display quality was not sacrificed to gain these advan- 
tages. Because of their shorter electron beam path length, 
magnetic-deflection CRTs give bright, crisply focused dis- 
plays. Vector scanning is limited in the number of traces 
that can be displayed before flickering becomes annoying. 
Raster scanning does not have this problem. This allows 
the HP 3561 A to display measurement results in three-di- 
mensional form with up to 60 traces of 400 points each 
(see Fig, 5 on page 10). That's the equivalent ol over 24,000 
vectors. Also, raster scanning consistently produces good- 
quality alphanumeric^. Well annotated displays <ire an im- 
portant part of the friendliness of the HP 3 561 A. 

Display quality is further enhanced by scanning the raster 
vertically at a rate of 36,205 Hz with the vertical lines 
formed from left to right. Since most rapid trace variations 
occur vertically, this gives a smoother trace by reducing 
the jagged line appearance compared to ei horizontally 
scanned raster. 

The digital portion of the display contains a memory 
control circuit, a special svn( circuit, and a processor inter- 
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face. The display is bit mapped with a vertical resolution 
of 256 pixels and a horizontal resolution of 512 pixels, The 
pixel bit map is stored in two sets of eight 1BK x 1 dynamic 
RAMs, One set is for full-intensity pixels and the other is 
for half-intensity pixels. Different Intensities are used to 
give emphasis to specific areas on the CRT. To allow the 
display to be synchronized to the power line frequency, a 
special sync circuit is used which does not require phase- 
locked loops or analog circuitry. The circuit free- runs dur- 
ing display updates and synchronizes during retrace inter- 
vals when minor timing variations are not noticeable. The 
addition of the synchronization circuit requires only one 
more flip-flop. 

The processor interface allows quick generation of dis- 
play graphics, It includes alphanumeric generation capabil- 
ity and simple line-drawing hardware that can draw verti- 
cal, horizontal, 45 c t and dotted lines at l r 400.000 pixels 
per second. It also supports several block operations such 
as set, clear, and complement at 11 ,300,000 pixels per second. 
These hardware features significantly reduce the amount of 
software overhead required for a bit -mapped display, 

Power Supply 

Requirements for the HP 3561 A's power supply were 
that it be lightweight and compact for inslru merit portabil- 
ity and that it provide eight different outputs, four of which 
are isolated floating supplies. A switching supply design 
seemed to be the natural choice for the job. The potential 
disadvantage of using this type of supply is that it can be 
a source of EMI [electromagnetic interference). For this 
reason, an out-of-band switching frequency of 128 kHz syn- 
chronized to the system clock is used. 

A half-bridge* off-line circuit topology was chosen, using 
500-vnlt VMOS switching transistors [Fig. 5). Control loop 
circuitry regulates the main 5V logic supply. The output 
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voltage is normally under control of a voltage feedback 
loop* but during overcurrent conditions a separate loop 
takes control to provide fold back current limiting. Both 
loops have their own loop shaping networks. All other 
supplies are derived from other windings on the same trans- 
former and l bus are slave preregulated by the regulation 
of the 5 V supply. Each of these other supplies is then further 
regulated by a linear three-terminal regulator. This pro- 
vides additional ripple rejection and supply isolation and 
the preregulation maximizes the efficiency of the linear 
regulators, 

The power supply is protected by four separate protec- 
tion circuits. The four fault-detection circuits are floating, 
nonfloating. and 5V overvoltage comparators and a pri- 
mary-overcurrent-sense circuit, All protection circuit out- 
puts are monitored by a latching shutdown circuit that turns 
off the gate drive to the powder VMOS transistors and indi- 
cates on one of four LEDs which fault condition occurred. 

Packaging and Shielding 

The HP 3561A package is similar to that used for portable 
oscilloscopes. The circuit boards sit in a vertical card cage 
and plug into a motherboard. Since the display and the 
switching power supply were identified as potential EMI 
sources, they are both placed behind a shield that runs the 
length of the instrument. 

Suppression of residual spurious signals by more than 
-80 dB relative to full scale required shielding of both the 
display and input amplifier. Electromagnetic field frequen- 
cies of concern are the display's vertical scan rate of 36 T 2Q5 
Hz and its harmonics. Electric fields caused by the high 
accelerating potentials for the CRT are attenuated by a 
grounded screen that covers the CRT face. High- impedance 
circuits in the HP 3561A"s input amplifier are covered by 
a grounded aluminum shield can. 
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Fig. 5, Block diagram of the HP 3561 A's power supply. 
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Varying magnetic fields are generated by the display's 
magnetic elements [flyback transformer and yoke) and by 
other ac high-current carrying conductors, These fields can 
generate currents in ihe instrument frame, causing local 
field radiation near the input amplifier. The initial field 
levels near the input amplifier were small (rniJHgauss 
levels), but the performance specifications required addi- 
tional attenuation of more than 40 dB. Shield materials 
have lower permeabilities at these low field levels, making 
the task more difficult. To simplify the problem, field radi- 
ation is reduced at the flyback transformer by a shorted 
turn of aluminum wrapped around the transformer, Cir- 
culating currents in this loop create opposing fields (Lenz's 
law), to help cancel some of the radiation, which in turn 
reduces field levels at the input amplifier. Further attenu- 
ation is achieved by enclosing the input amplifier board 
in a MuMetar* shield. The CRT assembly (including yoke) 
is enclosed in another MuMetal shield, both to reduce its 
radiation and to shield the CRT from external fields thai 
can cause display distortion. Fields from the power sup- 
ply's magnetic components have frequencies at 128 kHz 
and its harmonics and thus are above the HP 356lA's mea- 
surement baseband of 100 kHz (one of the primary reasons 
for switching the power supply at 128 kHz], Aliased spuri- 

MuMeEaf is a 'raoemarK of Allegheny- LlkHu«) Sieel Corporation. 



ous suppl}' radiation problems are eliminated by syn- 
chronizing the supply to the system clock, placing the 
aliased products at dc, Careful attention to grounding tech- 
niques also played an important part in the total solution. 
The resulting spurious signal levels are less than — 80 dB 
relative to full scale on the HP 3561 A 's most sensitive range 
51 dBVrms). 
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Instrument Software for Dynamic Signal 
Analysis 

by Glenn R. Engel and Donald R. Hitler 



THE SOFTWARE for a microprocessor-controlled in- 
strument determines the instrument's usefulness 
and personality — how it interacts with the user and 
how easy it is for a new user to understand and use effec- 
tively. Because there are many variables in dynamic signal 
analysis measurements, the software for the HP 3561 A 
Dynamic Signal Analyzer is particularly important. Some 
of the key features provided by the IIP 3 561 A "s operating 
system are: 

■ Friendly ustrr interface 

■ Autocall oration 

■ Overlapped processing for quick display updates 

■ Spectral map display 

■ Direct output to HP-IB (IEEE 488) printers and plotters 

■ Nonvolatile memory storage 

■ Self-tests and service tests. 

Architecture 

At the heart of the HP 3561A"s software is a 4 K* byte 
multitasking opera ling system that provides process syn- 



chronization and resource control. Thare nn live operating 
system tasks in the HP 3561 A (see Fig. 1). 

■ I/O control 

■ Command processing 

■ Measurement/display 

■ Nonvolatile memory con! ml 

■ Message handling, 

The I/O control process interprets commands from the 
HP-IB and front panel. This process checks syntax and 
parameter limits and does immediate execution of simple 
commands. An example of a simple command would be 
setting the marker position. Commands that cannot be exe- 
cuted immediately such as buffer dumps and nonvolatile 
memory operations are passed on to the command process 
for completion. The command process then coordinates 
the operations necessary to complete the command. This 
includes locking necessary resources t communicating with 
other processes lo handle parts of the command, and then 
unlocking resources after command completion. 

The measurement/display process is responsible for all 
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signal processing functions and display updates. Signal 
processing functions include data collection, fast Fourier 
transforms of data (see box on page 20], amplitude correc- 
tions, phase corrections, averaging, and data scaling. The 
measurement/display process takes advantage of over- 
lapped processing to increase data throughput. This pro- 
cess is shown in Fig. 2 and is discussed later. 

The nonvolatile memory process is responsible for all 
nonvolatile memory functions such as saving and recalling 
traces and instrument states. With the optional bubble 
memory installed, the nonvolatile memory process be- 
comes very important, because of the slower access time 
of bubble memory devices. By handling memory operations 
with a separate process, the HP 3561 A can be making useful 
measurements while traces or states are being saved or 
recalled from the bubble memory. 

The message process handles error messages and user 
warnings, This function Is a separate process because mes- 
sages need to be prioritized and most warnings have a 
display timeout. By using the message process, other pro- 
cesses do not have to worry about what happens to a mes- 
sage once the message process knows about it. 

The operating system provides several functions to ma- 
nipulate semaphores. The primary request and release func- 
tions provide a synchronization mechanism for resource 
allocation and are used to control counting semaphores for 
message communication between processes* When a pro- 
cess requires a resource 1 it calls the operating system's 
request function. If the resource is available, the process 
continues execution. If not, the operating system suspends 
the process and does not allow it to run until the resource 
is available, A resource becomes available when the process 
using it calls the operating system release function. When 
a resource is released, the operating system checks for any 
processes waiting for the released resource and schedules 
them for execution. 

Friendly User Interface 

The software for the HP 356 1 A 's user interface was writ- 
ten to provide the user with as much feedback as possible 
concerning error conditions and instrument operation. For 
example* if the user enters an illegal value for a numeric: 
input, the HP 3561 A, instead of just beeping, responds by 
telling the user exactly what the mini muni and maximum 
values can be. In addition, the software helps the user by 
giving suggested operation hints. For example, when the 
pseudorandom noise source is selected, the HP 3 5 61 A re- 
minds the user that the uniform window is the most likely 
window selection for that noise source. If the uniform win- 
dow is already chosen, no message is given. The friendly 
user interface extends to the HP -IB interface as well as 
front- pan el operation. A diagnostic mode a] lows for easy 
debugging of HP -IB error conditions. This mode displays 
each character received from the HP-IB interface on the 
instrument's display and marks the exact location of syntax 
violations, which can save programmers considerable de- 
bugging time, 

Autoealibration 

Calibration in the HP356lAisa three-step process, First, 
the software characterizes the amplitude and phase re- 



sponse of the instrument's analog front end over the entire 
operating range, Second, the software computes \)m front- 
end correction for the frequency range currently selected 
for measurements. This range, or span, can be a small por- 
tion of the HP 3561A's baseband of 100 kHz. Third, the 
software simulates the digital filters to determine the effects 
of the filtering operation on the magnitude and phase of 
the input signal. This results in one correction that repre- 
sents any errors introduced by the HP 3 561 A between its 
input connector and the final digital representation of the 
input signal in RAM. 

The key to providing high amplitude accuracy for the 
HP 3 5 bl A is directly related to its ability to measure its 
own magnitude and phase response properly. An accurate 
digital-slate-to-analog-level [1-bit D-to-A] converter deter- 
mines the accuracy of the internal calibration signal; it is 
driven by a pseudorandom noise generator. The spectrum 
nl j Ins calibration signal consists of a number of frequencies 
(determined by the length of the pseudorandom sequence) 
that are of nearly constant magnitude and whose magnitude 
and phase are theoretically exactly known. This allows 
rapid calibration cycles, since both magnitude and phase 
of all desired calibration frequencies can be determined 
with a single operation. 

One problem of this autoealibration scheme is slew-rate 
limiting of the analog version of the signal, since the ideal 
signal has step changes between the two levels. If the slew- 
rate limiting is symmetrical, the effect is that of passing 
the signal through a low- pass filter with an impulse re- 
sponse that has the width of the slew time. For high slew 
rates (which are present in the HP 35f$IA*s calibration cir- 
cuit], this filtering causes negligible errors, On the other 
hand, if the slew rates are not symmetrical, the effect is 
that of adding a small signal that is the derivative of the 
pseudorandom sequence. This additional signal causes 
small but significant errors in the amplitudes of the 
pseudorandom noise frequencies and the resulting calibra- 
tion constants. Fortunately, these errors are nearly identical 
and of the same phase when an amplitude-inverted 
pseudorandom sequence is passed through the same cir- 
cuitry. This allows the error to be canceled by subtracting 
a measurement of the inverted sequence from a measure- 
ment of the non inverted sequence. Another way of looking 
at this situation is that the errors are even-order and can 
be cancelled by the same technique used in a differential 
transistor pair or a push-pull amplifier to reduce even-order 
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Fig. 1. Operating tasks for the HP 3561 As software. 
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Fig, 2* (a) The result of serial pro- 
cessing— data is lost while the 
processor is computing the FFT 
and updating the display By 
using overlapped processing (b), 
records can be processed at a 
greater rate without losing data. 



distortion. Differences in propagation delays in the digital 
circuitry driving the analog conversion circuitry also cause 
the same kind of errors and so are also canceled. 

By using pseudorandom noise autncalibratiom the HP 
3561 A is able to correct for antialiasing filter rnlluff well 
past the cutoff frequency of 100 kHz. However, because of 
possible aliasing problems and decreased dynamic range, 
the accuracy tif the HP 3561 A is not specified above 100 
kHz. The resulting accuracy specification is ±0.15 dB from 
to 100 kHz. 

Overlapped Processing 

Overlapped processing is possible in the HP 3561 A 
thanks to a multipart RAM. This RAM is accessible from 
four different sources: the main processor, the slave FFT 
processor, and two DMA [direct memory access) channels. 
This allows Ihe HP 3561 A to almost triple its real-time rate 
by dot Jig ihree operations simultaneously. Fig. 2 illustrates 
the possibilities. While the main processor is manipulating 
record N for display, the FFT processor is processing record 
N + l, and the DMA sequence la collecting record N + 2. 

Spectral Map Display 

The software for the HP SSGlA's spectral map display 
[see Fig. 5 on page 10) is quite simple. This display format 
provides a three-dimensional representation of the analyzed 
input signal versus time. As new spectra are placed into 
the map, they must be drawn as if Ibcv were behind I he 
previous spectra. This operation is similar to viewing a 
distant range of mountains from an airplane. Mountain 
ranges that are behind others are only visible if then are 
higher than the ranges in front of them. That is, higher 
mountains can obscure lower mountains he bind them. The 
hidden-line algorithm used by the HP 3561 A uses exactly 
this idea. When a new spectrum is put on the display it is 
compared point by point with the previous spectrum to 
decide whether it is ** visible/' This process is simplified 



for a raster scan display because the two spectra only need 
to be compared at a limited number of points — once every 
horizontal pixel position. No interpolation is needed be- 
tween points and only the largest value for each horizontal 
pixel position DMisi be remembered for later spectra, 

HP-IB Plots and Printer Dumps 

HP- IB f>lnts and printer dumps are supported on the HP 
:J5ti 1 A to allow the user to produce a hard-copy represen- 
tation of the display. Printer dumps are achieved by simply 
reading the bit-mapped display memory and sending this 
information to a printer that supports raster-type graphics. 

Plots are obtained by sending HP-GL plotter commands 
to a pen plotter. A high-quality copy of the HP 3561 A 
display is obtained in this manner. Plotting the spectral 
map is more complex than generating the raster-display 
spectral map, When plotting the spectral map, the HP 
3561A must keep a copy in its memory of the largest values 
plotted to aid in hidden-line removal. These values are 
kept with enough resolution to provide correct hidden-line 
removal down to the width of the plotting pen. HP-IB plots 
also allow the user to annotate specific marker values to 
show test results or merely to provide extra documentation* 
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FFT Implementation 



One of the benefits of using an analyzer based on a fast Fourier 

transform (FFT) is measurement speed. The HF 3561 A Dynamic 
Signal Analyzer uses the FFT algorithm to analyze signals with 
the equivalent of a parallel bank of 400 narrowband filters. This 
simultaneous spectral analysis technique is inherently much fas- 
ter than using a swept-trequency spectrum analyzer that must 
cover the same frequency span serially, one frequency at a time. 
Although not yet practical at RFand microwave frequencies be- 
cause of the fast FFT computation rate required, current technol- 
ogy does aifow its use at low frequencies This speed improve- 
ment can be put to good use in many of the HP 3561 As appli- 
cations The typical display update rate of five to ten spectra per 
second makes circuit adjustments easier A nojsy signal can be 
averaged 1000 times in less than a minute in the fast averaging 
mode. 

These fast measurements require a fast FFT computation. En 
the HP 3561 A. the FFT computation is performed by a TMS320 
microprocessor. This processor has a reduced instruction set 
designed for digital signal processing applications The proces- 
sor and its support chips occupy only 30 square inches of printed 
circuit board area. It computes a 512-poinl complex FFT in 35 
ms and a 1024-point real FFT in 42 ms. including the window 
multiplications. 

Speed 

Fast execution times are obtained by tailoring the FFT algorithm 
to the internal architecture of the TMS320. A modified radix-8, 
Sande-Tukey twtddte-factor algorithm 1 is used. With a slight re- 
arrangement of the multiplications, this algorithm fits nicely mto 
the TM5320"s sum-of-producls architecture. A radix-8 algorithm 
was chosen for several reasons With this radix, bit-reversing is 
replaced by digit-reversing, that is, groups of bits are reversed 
at one time, which is easier and only needs to be done for every 
eighth data address. Also, a radix-8 algorithm requires iewer of 
the slower I/O operations to the instrument's main RAM, To reduce 
the effect of slow RAM accesses further, the algorithm performs 
arithmetic operations using the on-chip cache RAM while waiting 
for the I/O operations 

Noise and Accuracy 

Noise performance of the FFT algorithm is equally important 
for the HP 3561 A, A numerical aigorithm is only as accurate as 
the number of bits used in the calculations Inaccuracies in cal- 
culating the FFT can degrade absolute amplitude accuracy and 
linearity and generate phantom spectral components These are 
all lumped into the category of "noise." In the MP 3561 A, the FFT 
noise performance can be no better than the RAM data bus 
width, which is 16 bits. Given a perfect algorithm (no errors), the 
maximum dynamic range would be approximately 96 dB (6 dB 
per bit). This is more than adequate to ensure the instrument 
dynamic range specification of 80 dB. including all scaling over- 
heads and specification margins, 

The challenge js keeping a 16-bit FFT aigorithm accurate to 
16 bits at its output The algorithm reads data words from main 
RAM. performs arithmetic operations, rounds the results and 
stores them back in mam RAM. This occurs once for each FFT 
level, and there are four or five levels per transform If a 16- bit 
data point is multiplied by a 16-bit coefficient, the result is 31 
bits long If this result is then summed with eight others just like 
it, three more bits are added and the result could overflow the 
TMS320 s 32-bit accumulator Therefore, prescaling is necessary, 

The HP 3561A's FFT algorithm uses a simplified block scaling 
technique It <s impfemented by saving the magnitude of the 



largest data point observed for an FFT level The next level 
chooses to divide by one, two. four, or eight based on the previous 
level's maximum. This scaling method gains much of the perfor- 
mance improvement of a true block floating-point algorithm with- 
out the added complexity or speed reductions. With this type of 
scaling, the instrument nosse floor drops with the level of the 
input signal. 

Another important consideration is the nature of the noise floor. 
A noise floor that is uniform versus frequency, or white, is pre- 
ferred. This means thai the average power of the error in each 
of the output frequency bins is the same. Most of the bins are 
the result of successive multiply-accumulate operations, where 
many bits are discarded in the rounding operation. It is fairly 
easy to keep the error in these bins consistent. The problem bins 
are the ones that are the result of an accumulate only In these 
bins, very few bits are thrown away during rounding. The error 
has only a few possible values, and there is an expected error, 
or bias, in the result. For example, if two 16-bit values are added 
and the 17-btt result is rounded back to 16 bits, the error is either 
zero or y s LSB (least-significant bit). There is a bias of % LSB. 
When these biased results are summed in following levels, the 
bias error adds up in a linear fashion. 

Generally. FFT noise floors have spur-like bumps located on 
the output bins with large expected errors A user might be 
tempted to interpret these as spurious components in the input 
signal even though they are outside of the specified dynamic 
range. To solve this problem, two types of rounding offsets are 
used. One results in a positive bias, and the other results in a 
negative b*as. By appropriately choosing between these when 
rounding, the bias errors cancel 

The FFT limits amplitude accuracy in yet another way: the 
windowing operation. The input data block is multiplied by a 
tapering function, or window, that forces the start and The end 
of the record to zero. This reduces the effects of leakage, which 
occurs when the signal is not periodic in the measurement inter- 
val, Windowing is often thought of in the frequency domain as 
the frequency response of the analysis filters. To make accurate 
amplitude measurements, the effective filter frequency response 
should be very flat over a frequency span of plus or minus half 
of a bin spacing The HP 3561 As flattop window function is flat 
to ±0.003 dB, much flatter than those used m other FFT-based 
analyzers. This flatter window is needed to complement the im- 
proved amplitude measurement accuracy of the HP 3561 A, 

Circuitry 

The TMS320 hardware looks like a slave processor to the HP 
3561 As 68000 main system processor The 68000 sets up a 

small table in main RAM that tells the TMS320 what type of FFT 
to perform, scaling restrictions, data buffer start addresses and 
other setup information. Then the TMS320's reset line is set true 
and FFT computation begins. When finished, The TMS320 writes 
the scale factors in main RAM and interrupts the 68000. Main 
RAM is shared among several processors and DMA ports. A bus 
arbitrator establishes priority and grants RAM accesses. 
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Kolesar 

A Compact, Reliable Power Supply for an Advanced Desktop Com- 
puter, Jack L, Burkman, Howell R FelsenthaJ, thomm 0. 
Meye r , an d Wt i rre n C P m 1 1 

Compact 32-Bit System Processing Units* Kevin W. Allen, Paul 
G Christo/anelii, Robert E, Kuseski, Ronald D. Larson, David 
Mai t land, and Larry J. Thayer 

June 1984 

A Parametric Test System for Accurate Measurement of Wafer- 

Stage ICs, Yoh Narimotsu and Keiki Kana/uji 
Powerful Test System Software Provides Extensive Parametric: 

i rumen t Capability, Takuo Banno 
A High-Speed 1 -Ml {z Capacitance /Conductance Meter for Measur- 
ing Semiconductor Parameters, Tomoyuki Akiyama and Kenzo 

Jshiguro 
An Electronic Tool for Analyzing Software Performance, Gail E. 

Hamilton, Andrew J. Blasciak, Joseph A. Hawk, and Brett K, 

Carver 
Counter Module Simplifies Measurements on Complex Wave* 

forms. Donald J. Smith, Johnnie L Hancock, and Thomas K, 

R oh ley 

July 1984 

A New Handheld Computer for Technical Professionals, Susan L, 

U'echsler 
Soft Configuration Enhances Flexibility of Handheld Computer 

Memory, Nathan Meyers 
Custom CMOS Architecture for a Handheld Computer. James P. 

Dickie 
Packaging the HP-71B Handheld Computer. Thomas B, Lindben> 
Module Adds Curvo-FHtiujj tind Optimization Routines to theH7lB, 

Stanley M. Rlascotv, jr. and James A. Donnelly 
ROM Extends Numerical Function Set of Handheld Computer, 
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Laurence W. Grodd and Charles M. Pa t ton 
Plug-In Module Adds FORTH Language and Assembler to a Hand- 
held Computer. Robert M. Mi/ler 

Augusl 19B4 

Touchscreen ParsoneJ Clomp titer Offers Ease of Use and Flexibility, 

Srlnivos Sukumor 
Operating System and Firmware of the HP 150 Personal Computer, 

Laurie E. Foliero Wood and Charles H. WheJan 
1 be ]IP 150 TcHji Iim rvrn. Aei I nl tractive User Input Device for a 

Personal Computer, Peter ll. Straton, Scott R. McClelland, and 

Thomas £\ Kiibourn 
Applications Software for the Toui Persona! Computer, 

PeterS, Showman, Kar J IV. P^tfjs, Kuril v / Arkm. Jeffrey A. Spoel- 

stm, John Price, W* Bruce Cufberlson, and Robert D. ShurtJeff. Jr. 
Hardware Design of tbe HP 150 Personal Computer, John E. Wctk j i i >■. 

Patricia A, Brown, George Szemnn, and Susan E. Carrie 
Software Graphics in (he HP nu 

Personal Computer Printer Is User Installable, Joseph D. Barbara 
A Standard Keyboard Family for HP Computer Products, Lorenzo 

Dunn und Michael FL Perkins 



September 1984 

Transmission Impairment Measuring Set Simplifies Testing of 

Complex Voice and Data (Jin .nits. David R. Novotny, Jeffrey 

Tomberlin* Charles P. Hill, fames P, Quan, Gordon A, lensen, 

and Jerry D, Morris 
TIMS Mechanical Design 

Weight, Size, and Noise Impact Power Supply and Display Design 
Master 'Slave TIMS Operation Improves Productivity, Teresa L, Reh 
How Master/Slave Mo do Works 
Testing the TIMs, Allan IV. Dodge, Scott S. NeaJ f and Kurt R. 

Goldsmith 
Semiconductor Research Corporation: A Perspective on Coopera- 

live Research, Richard A, Lucie 
A Hyphenation Algorithm for HPWord. PtiuJ R. Smil 
Designing Software for the International Market, Heather Wilson 

and Michael J. Shaw 



October 1984 

The HP 3065 Board Test Family: A System Overview, 

R. Fay and John E. McDermid 
HP Q-STAR 



Thomas 



Confirmation-Diagnostics 

Automatic Test Program Generation for Digital Board Testing, 
Robert E. Boilievv 

Board Test Connection Terminology 

Digital Subsystem for a Board Test System, Matthew L. Snook 
and Michael A, Tesrka 

Digital Test Throughput 

Safeguarding Devices Against Stress Caused by In-Circuit Test- 
ing, Vance R. Norwood 

Extensive Library Simplifies Digital Board Test Setup, Bandy W. 
Ihilmherg 

An Interpreter- Based Board Tesl Programming Environment, 
Mark A, Mulhirw 

Testing for Short-Ciriuit Failures. T Michael Hendricks 

Reducing Errors in Automated Analog In-Circuil Test Genera- 
Hon, John E. McDermid 

November 19B4 

An Advanced 5-Hz-to- 200-MHz Network Analyzer, Robert A. 

Witfe and Jerry W. Daniels 
User-Defined Vector Math Expand* Measuremenl Capabilities 
A Broadband Two-Port S-Pararaetei Tbb\ Set, William M, Spauiding 
An ADC for a Network Analyzer Receiver, Alan J. Baker 
An Industrial Workstation Terminal for Harsh Environments, Jean 

Bounaix, jean-Claude Dureau. and Jacques Firdmann 
How Do You Describe Terminal Ruggt:dness? 
High-Quality r Dot-Mai rixlniptJtJ Printer Family. Mark}. Di Vitrorio 
Custom IC Controls Dot-Matrix Printers, Thomas B. Prilchord and 

David S. Lee. 

December 1984 

Versatile instrument Simplifies Dynamic Signal Analysis at 

Low Frequencies, James S. Epstein 
Dynamic Signal Analysis for Machinery Maintenance 
Hard wan j DE^si^n for a Dynamic Signal Analyzer, James S. Epstein, 

Glenn R. EngeJ, Donald R. Hiller, Glen L. Purdy, Jr., Bryan C. 

Hoog, and Eric J. VVicklund 
Instrument Software for Dynamic Signal Analysis, Glenn R, Engel 

and Donald R. Hiller 
FFT Implementation 
Custom Digital Filters for Dynamic Signal Analysis. Charles R. 

Pa nek and Steven F. Kator 



PART 2: Subject Index 



Subject 

A 

Acoustic measurements 
Active input and output 



Month 



Dec, 



circuits, TIMS 


Sept, 


Actuator, disc drive 


Jan, 


ADC. dynamic signal analyzer 


Dec. 


ADC. receiver, network analyzer 


Nov. 


ADC. 16-bit. low-cost 


Apr, 


Algorithm, board test connections 


Oct, 


Algorithm, Booth 


Feb. 


Algorithm, Fh>tc;her-Powell 


Inly 


Algorithm, hyphenation 


Sef.it. 


Algorithm, short circuit testing 


Oct. 


Algorithm, touchscreen scanning 


Aug, 


Alignment, automatic. 




disc drive Etc 


fan. 


Amino acid analysis 


Apr. 


Amplifier, prescaling, ADC 


Nov. 


Analog testing, in-circuit 


Oct. 



Analyzer, dynamic signal 
Analyzer, logic timing 
Analyzer, oetwoxk 
Analyzer, software performance 
Architecture, 32-bil processing units 
Architecture, handheld computer 
Architecture, parametric test system 
Array red! mens ioning 
Assembly language ROM 
Aiil>H.,.,lihrriliim. dynamic: signal 

analyzer 
Autoinjector, LC 
Autos am pier, LC 
Averaging, time interval 

B 

Baluns, broad frequency range 
BASIC" programmable computer 
BASIC. 3-D graphics 
BASIC. Board Test 



Det: 


BASIC. PIP 9000 Model 520 


Mar, 


Feb. 


BASIC. I'O, HP 9000 Model 520 


May 


Nov. 


BASIC, multiprogramming 


May 


June 


Bed -of- nails lest fixture 


Oct, 


May 


Bell-standard transmission 




July 


measurements 


Sept. 


June 


BIOS 


Aug, 


luly 


Bit-slice data acquisition processor 


Apr, 


luly 


Board testing, in -circuit 


Oct. 




Bond-wire heat i tig, in-circuit testing 


Oct, 


Dec. 


Booster pump 


Apr. 


Apr. 


Booth algorithm 


Feb, 


Apr, 


Bridge, reflecto meter, configurations 


Nov, 


j Li lie 


Buliulejei 


Mar. 




Buffer, circular, IO 


May 



[uk- 
May 
Oct. 



Cable compensation* capacitance 

measurements June 

Cable conned itms. sealed Nov. 
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Cache memory Feb, 

Cache memory, 1 O Mar, 

CALC mode. HP-71B July 
Capacitance measurement 

subsystem June 
Card file software, personal 



computer 


Aug. 


Card reader/recorder, handpulled 


July 


Case. sealed 


Nov, 


Ceramic disc, rotary valve 


Apr. 


Clock recovery, disc drive 


Jan. 


CMOS, handheld computer system 


July 


CMOS testing, latchup prevention 


Oct. 


Code and data separation 


Feb. 


Code segment mapping 


Feb. 


Code, intermediate 


Mar. 


Codes, disc drive 


Jan, 


Coding, linear predictive 


Jan. 


Coherence breaking by phase 




modulation 


June 


Columns, liquid chromatography 


Apr. 


Command Set 8fl (for disc drives! 


Jan. 


Communications processor [COM) 


Apr. 


Communications units 


Apr, 


Compiler, runtime, interactive 


Mar. 


Complex data type 


Julv 


CompU.x mat hematics 


July 


Computer. 32-bit 


Max. 




May 


Computer, handheld 


July 


Computer, personal, LC controller 


Apr. 


Computer, touchscreen, persomil 


Aug. 


Computers, IIP lUUO A-Series 


Feb. 


Configuration, automatic 


June 


Confirmation/Diagnostics 


« >, i 


Connection modes, capacitance 




measurements 


June 


Connector, t ustnnj pnwor transistor 


May 


Contexts, SI IN operating system 


Mar. 


Control system, LC 


Apr. 


Cooling, desktop workstation 


May 


Cooling, power supply 


May 


Cooperative research 


Sept. 


Counter, gated universal 


lum: 


i-inJHj.iiun molding 


SepL 


,1 quartz, measurements 


Nov. 


1 .nrvn lilting 


Inly 


Custom IC, printer control 


Nov. 


Custom ICs. counter 


June 


Custom ICs. digital filter 


Dec. 


C-HIGH RESOLN function 


June 


C-1 measurement 


June 


C-V measurements 


June 


D 

Daisy-chain cmifigura! inn 


July 


Damping, LC 


Apr. 


Data acquisition processor (DAP) 


Apr. 


Data circuit testing, Bell-standard 


Sept. 


Data communications, workstation 


Mar, 


[Jdtti storage and retrieval, logii 




analysis 


Feb. 


Datacom drivers 


Aug, 


Debugger, software 


Mar. 


Degradation table 


SepL 


Detector, touchscreen 


Aug. 


Detectors, liquid chromatography 


Apr. 



Development, parallel 

Dibit integrator 

Diagnostic control system 

Diagnostics, 3 2 -bit SPUs 

Digital filter 

Digital subsystem, board lest 

Digital testing, in-circuit 

Diode array detector 

Disc drives 4 04 M -byte 

Dispersion, holographic grating 

Dispersion. LC 

Display, color graphics 

Display, dynamic signal analyzer 

Display, tilt mechanism 

Dithering, dynamic signal analysis 

Dot -matrix printer family 

Duty cycle measurements 



EMI design, rope gasket 
Environmental IP rating 
Environments, industrial. 

terminal for 
Error correcting memory 
Error correction, capacitance 

measurements 
Error detection and correction , 

disc drive 
Errors, analog in-circuit testing 
Events measurements 
Exchanges 



Fast Fourier transform 

Fasteners, tolerance design 
FFT analysis, low-frequency 
Pile management system 
I ilir manager, personal compute 

security 
Filter, digital 

Filterphotometric detector 
Financial calculator software, 
personal computer 

Fletcher-Powell method 

J l( i.iling-potnt chip set 

Flow cell 

FM code, Manchestei 

FORTH language ROM 

FORTH/BASIC conversion 

Fraclional'N synthesizer board 

l-'nu Mnruil-N synthesizer, modified 

inncy measurements 
Front panel, virtual 

G 

Gated universal counter 

Generic research 

Graphics software, personal 

computer 
Graphics. 3-D, BASIC 
Graphics, asynchronous input 

phics, input device tracking 
Graphics, text 

Grating) holographic diffraction 
Group delay, measuremiHil 



Mar. 


H 




Jan. 


Head alignment, automatic, disc drive Jan, 


Feb, 


Head positioning system, disc drive 


fan. 


May 


High-pressure pump 


Apr. 


Dec. 


High-speed columns, LC 


Apr. 


Oct. 


Holographic diffmction grating 


Apr 


Oct. 


Homer's method, polynomial roots 


July 


Apr 


HPLC 


Apr. 


Jan. 


HPYVord hyphenation algorithm 


Sept. 


Apr. 


HiMB transformer 


Apr. 




HP-1L, interface module 


July 


May 


HP Q-STAR, network, software 


Oct. 


Dec. 


HP-UX operating system 


Mar* 


May 


HP-UX. corporate strategy 


Mar 


Dec, 


Hyphenation algorithm 


Sept. 


June 


I 






IEEE floating-point standard 






In-circuit device testing 


Oct. 


May- 


Industrial workstation terminal 


Nov. 


No v. 


Infrared array touchscreen 


Aug. 




Injector, LC 


Apr, 


Nov. 


Tnk jet printing, thermal 


Mar- 


Feb. 


Integral function 


July 




Integrals, nesting 


luiy 


lum-i 


Interfaces. I/O, TIMS 

i Hi rrmod illation distortion, 


Sept. 


Jan. 


measurement circuit 


Sept. 


Oct. 


Interrupt handling 


M.-ir. 


June 


IF rating 


Nov. 


Aug. 


1PG-1I, automatic test program 






generator 


Oct. 




ISO substitution tables 


Sept. 


July 


[ tu-V converter 


lunv 


IM:. 


M) drivers 


Mar. 


May 


I/O subsystem, 32-bll computer 


Mai 


Dec. 


I/O. BASIC, examples 


May 


July 


I/O. TRANSFER statement 


May 


Aug. 


I/O, unified 


May 



July 
Mi-. . 
Apr. 

Aug. 

Jan. 

]nlv 

Feb. 

Apr. 

Jan. 

July 

July 

Nov. 

Sept, 

Jurm 

June 



lllllr- 

Sept. 

Aug. 
May 
May 
May 
May 
A pi 

Nov, 



Junction temperature, in-circuit 
testing Och 

K 

Keep/toggle vectors Oct 

Ki , redefinition Inly 

Keyboard family, computer products Aug. 
Keyboard, detached design 

Keyboard, sealed Nov. 



Laguerre's method, polynomial rools July 

LAN 9000. local area networking Mar. 

Level measurements. TIMS Sept. 

Library, digital tesl Oct, 

Library, speech Jan. 

Library. Waveform measurement Juup 
Lightweight frame, TIMS 

Limit cycles, digital filter Dte 

Linguistic rules Sept. 

Liquid chromatograph Apr. 

LO board, network analyzer Nov. 

Localization, software Sept. 

Logic development! timing analyzer Feb. 
Longitudinal balance measurement, 

TIMS Sept. 
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Low-dispersion liquid 



chro ma to gr a p h y 


Apr 


Low-pressure compliance 


Apr, 


Low-pressure pump 


Apr. 


M 




Machinery maintenance 


Dec. 


Man agem en t , project 


Mar, 


Manchester code, I VI 


Jan, 


Marked events, logic analysis 


Feb. 


Mrisk design, filter IC 


Dec, 


Master/slave mode. TIMS 


Sept. 


Matrix calculations 


July 


Measurement, dc subsystem 


June 


Mechanical design, :i 2-bit SPUs 


May 


Mechanical design, disc drives 


Jan, 


Mechanical design, power supply 


May 


Mechanical design, TIMS 


Sept. 


MfirioMaker 


Aug. 


Memory, overlay structure 


Feb, 


Memory, soft -configuration 


July 


Memory, virtual 


Mar. 


Mentors , industrial 


Sept, 


Metering pumps 


Apr. 


MFM code 


Jan. 


Microbnre columns, LC 


Apr. 


Micromachine, HP 100U A-Series 


Feb, 


Microparaphraser 


Feb. 


Mixer-multiplier, digital filter 


Dec. 


Mixing chamber 


Apr. 


Mobile phase 


Apr. 


Modeling, foam core board 


May 


Modular design, desktop workstation 


May 


MS-DOS operating system 


Any. 


Multiple-processor system design 


Mar. 


Multipoint conrmc(icms< industrial 




terminal 


Nov. 


Multiprogramming, BASIC 


May 


Multiuser facility 


Feb. 


Mute character table 


&Bfft 


N 




Narrowband measurements , 




low- frequency 


Dec, 


Network analyzer 


Nov. 


Networking, board tesl 


Oct. 


Networking, local area 


Mar. 


NLQ h near letter quality printing 


Nov. 


NMOS-lll ic;s 


May 


Noise reduction, TIMS display and 




power supply 


Sept. 







Octave analysis, acoustic 




measurements 


Dec. 


Offset errors, digital filter 


Dec, 


Offset frequency board 


Nov, 


Opens test 


Oct. 


Operating system, dynamic signal 




analyzer 


Dec. 


Operating system, HP-UX 


Mar. 


Operating system, MS-DOS 


Aug, 


Operating system , RTE-A 


Feb. 


Operating system, SUN 


Mar. 


Optimisation, example 


July 


Oscillator, 100-MHz time base. 





phase modulaled Jumh 

Oscilloscope measurement system June 

Oven, crystal reference Nov. 

Oven, LC Apr. 

Overlapped processing Dec. 



Packaging, handheld computer July 

Parametric testing, semiconductor June 

Peak-to-average ratio measurements Sept, 

Period measurements June 

Personal Applications Manager Aug, 

Persona] computer, touchscreen Aug. 

Phantom short circuits, testing Oct. 

Phase modulation, random oscillator June 

Ph.-tM' shift measurements June 

Phonetics, speech synthesis Jan. 
Pipeline shift registers. 

digital f 1 1 1 r ■ r Dec. 

Pipelined data path Feb. 
Point-to-point connections, 

industrial terminal Nov, 

Polynomial root finder July 

Postprocessing, logic analysts data Ki-k 

Power supply design, dual Sept. 

Power supply, desktop computer May 
Power supply; dynamic signal 

analyzer Dec 

Precision machining Apr. 

Prehealer, LC mobile phase Apr. 

Printer, t hermal , personal computer Aug. 

Printers, dot matrix Nov. 

Print head, carriage control Nov. 

Print wires, drive control Nov. 

Prior invention, coping with Mar. 

Probing pattern generator June 

Processing unils. 32-bit system May 

Production tests, TIMS Sept. 

Prog j 'a mining, board test Oct 

Propagation delay measurements June 

PTH amino acid analysis Apr. 

Pulse widtli measurements June 

Pumps, metering Apr. 

Pump, booster Apr. 

Purity, LC peak Apr. 



RACE program Sept. 

Rating, IP Nov. 

Read write electronics, disc drive Jan. 

Real-time computers, 1-3 MIPS Feb. 

Receiver, digital, IF Nov. 

Receiver, tims Sept. 

Reference board, frequency Nov. 

Regulator, primary switching May 

Reliability assessment tests, TIMS Sept, 

Removable disc pack, 4 04 M -byte ]iti. 
Research, semiconductor, genpri< . 

cooperative Sepl. 

Rise/fall time measurements hine 

Rotary valve. LC Apr. 

RTE -Aoperatingsystem Peb. 

Rugged industrial terminal Nov. 

Run-time compilation Mar, 



Safeguards, in-circuit lest OcL 



Sampler, LC Apr. 

Sampling, software measurements June 

S-paiameter measurements Nov, 

Sapphire metering pump pistons Apr. 

Scanned infrared array Aug, 

Scanner card, digital Oel. 

Sealing, industrial terminal Nov. 

Sector sparing Jan. 
Semiconductor parameter 

lesting June 
Semiconductor Research 

Corporation Sept. 

Separation, code and data Feb. 

Shorts lest i h I 
Signal analysis, dynamic, 

low-frequency Dec. 
Signed digit encoding Feb. 
Softkey menus, network analyzer Nov. 
Software graphics Aug. 
Software localization (translation] Sept. 
Software translation Mar. 
Software. C/C-V Meter June 
Software, digital analyzer Dec. 
Software, multiple-processor system Mar. 
Software, parametric test system June 
Software, performance analysis June 
Software, preserving investment Mar, 
Software, speech output Jan. 
Software ► touchscreen computer Aug, 
Software , waveform measurement June 
SOLVE function July 
Solvent delivery system, LCI Apr- 
Spectra, LC Apr, 
Spectral map display Dec. 
S pec trophoto metric detector Apr. 
Speech output, HP 1000 Computer Jan. 
Speech output, HP 3000 Computer Jan, 
Speech output, Series 80 Computers Jan. 
Spreadsheet program, personal 

computer Aug. 

Spring, micro Apr, 

Stationary phase Apr, 

Statistics, software analyses June 

Statistics, time interval Feb, 

Strip-line design, broadband Nov. 

Synthesis, frequency Sept. 

Nov. 

Synthesis, speech Jan. 

Stimulus/measurement units June 

SLIN operating system Mar. 

Supervisor board, power supply May 

Switching matrix, lest head June 

Switching, "dry' 1 technique June 
Symbolic interface, software 

analyzer June 



Tasks Aug. 

Temperature control. LC column Apr, 

Terminal nppratin^ system (TQS) Aug. 

Terminal, industrial Nov. 

Terminals, character set differences Sept. 
Terminology, board test 

connections Dei. 
Test instruction set, 40 62 A Test 

System June 

Test station, power supply May 
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Test system, in-circuit board 


Oct. 


Touchscreen personal computer 


Aug, 


Test system, parametric 


June 


Track follower, disc drive 


Jan, 


ng> automatic program 




Transactions, disc drives 


fan. 


generation 


Oct 


Transient analysts 


Dec 


Testing. EMI. computer 


May 


Translation, software 


Sept. 


Test in g, en vi ra nmen * a L computer 


May 


Transmission Impairment 




g, environmental. TIMS 


Sept. 


'j ring Set 


Sept. 


-n-circuit 


Oct. 


Transmitter. TIMS 


Sept. 


Testing, safeguards 


Oct. 


Tree, linguistic rules 


Sept 


Testing, short-circuit 


Oct. 


Triggers, counter, main and 




Test programming, interpreter-based 


Oct. 


delayed 


June 


Thermal printer, personal computer 


Aug, 


Two-counter time interval 




Thinkjei 


Mar. 


measurement 


June 


Throughput, digital testing 


Oct. 


Twisted-pair test wiring 


Oct. 


Time capture measurements 


Dec. 






Time interval averaging 


June 


u 




Time interval measurement 


June 


Universal counter gated 


[une 


Timing diagram, logic analysis 


Feb, 


UNIX, implementation. HP 9000 




TIMS 


Sept. 


Computers 


Mar, 


Iili-ranced design 


\lav 


User microprogramming 


Feb, 



UV/Vis detectors. LC 



Apr- 



Valve, rotary, LC 


Apr, 


VCL. vector control language 


Oct. 


Vector math, network analysis 


Nov 


Vibration measurements 


Dec. 


Virtual control panel 


Feb. 


Virtual front panel 


June 


VLFM code 


Jan. 


Voice circuit testing. Bell-standard 


Sept. 


W 




Waveform measurement library 


June 


Wire frame. TIMS 


Sept. 


Word processing software 


Aug. 


Workstation* engineering 


Mar. 




May 



Zoom operation, digital filter 



Dec. 



PART 3: Model Number Index 



HP-71B 


Handheld Computer 


July 


HP 9000 


HP Series 80 


Computers 


|hii. 








Apr. 


HP 9000 


HP 150 


Personal Computer 


Aug. 




A600 


Computers 


Feb, 


HP 9000 


A700 


Computers 


Feb, 




A900 


Computers 


Feb. 


16320A 


HPIOOOA-Senes 


Computers 


Feb, 


igar.H.iA a 


HP 1000 


Computer 


Jan. 


19801 A/8 






Oct. 


2720 I A 


1040 A 


HPl ' ii System 


Apr, 


272(1 :.\ 


1090 


Liquid Chromatograph 


Apr, 


27205A 


134 5 A 


Digital Graphics Display 


Nov, 


2934DS 


1965 A 


Gated Univer$i4 Count© 


lnm- 


;{T 3 fi77A/B 


1980A/B 


i )m illoscope Measurement 




4 5405 A 




System 


June 


45410A 


2199C/D.2439.2139A 


A 900 Computer 


Feb. 


45420A 


^1!J7C/D P 2437.2137A 


A 700 Computer 




45422A 


2107 AK 




Feb. 


45423A 


21 96013,2436, 2 136 A 


A600 Computer 




46010 


2 100 A K 




Feb, 


64000 


2674A 


Thermal Printer 


Aug. 




2032A 


General Purpose Printer 


Nov. 


643 10 A 


2933A 


Factory Data Printer 


Nov. 


B4fiOOS 


2934A 


Office Printer 


Nov. 


7<W35A 


HP 3000 


Computer 


Jan. 


79846A 


306fi 


Hoard Test Svstem 


Oct. 


79847A 


3065C 


Controller 


On i 


79880A 


3O05H 


Test Station 


Oct. 


79881 A 


3081A 


J nd ustrial Workstation Terminal 


Nov. 


H2400A 


3501 A 


Dynamic Signal Analyzer 


Dec. 


a^401A 


3577A 


Network Analyzer 


Nov. 


82420A 


4062A 


Semiconductor Parametric Test 




82441 A 




System 


fune 


824H0A 


4084A 


Switching Matrix Controller 


June 


82484A 


4085A 


Switching Matrix 


June 


82967A 


4141A 


DC Source/Monitor 


June 




4280A 


l-MHzCMHt^C-VPloM<M- 


June 


9218B 


4945A 


Transmisskra Impairment 




97050A 




Measuring Set 


Sept. 




7933 


Disc Drive 


Jan. 


97052A 


7935 


Disc Drive 


Jan. 




HP 9000 


Sri Uts2(l(J( ;.i3iifmrrrs 


Mar. 


97070A 


HP'inmi 


Series 500 Computers 


Mar. 








May 


97079A 


HP 9000 


Model 520 Qimpi Mi i 


May 





Model 530 System Processing 

Unit May 
Model 540 System Processing 

Unit May 
Detached Keyboard, 

HP 9000 Model 520 Computer May 

Pin Board June 
Waveform Measurement Library }une 
Waveform Measurement Library June 

Speech Output Module Jan, 

Speech Library Jam 

Speech Library Ian. 

Single Bin Sheet Feeder Nov, 

S - P arame te r Tea t Sets Nov . 

VisiCalc*150 Aug. 

Series 100/Graphics Aug. 

MemoMaker Aug. 

Personal Card File Aug. 

F i n anc i a I Calcu I at or A ug. 

K*;y board Aug. 

Logic Devel opment System Feb. 

June 

S i if l Wrtre Perform a nee Analyzer (line 

Logic Timing Analyzer Feb. 

Solvent Delivery System Apr. 

Autoinh ;<[ur Apr. 

AutosampleT Apr, 

HPLG Detection System Apr, 

Filterphotometric Detector Apr. 

Card Reader Module July 

HP-IL Interface Module July 

4 K M e in o r y M od u 1 e July 

FORTI L Assembler Pac July 

MathPac July 

Curve Fitting Pac July 
Speech Synthesis Module, 

Series 8€ Jan, 

CbaracterFonts Nov, 
BASIC Language System. 

HP 9000 Model !i 20 Mar. 
3D Graphics, HP 9000 

Model 520 May 
HP-UX, HP 9000 

Model 520, single-user Mar. 
HP-UX , HP 9000 . Models 

530 and 540, single-user Mar. 
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07oauA 


HP-UX. HP 9U0Q 




9B760A 




Model 520 , mu 1 tiuser 


Mar* 


BH770A 


97089A 
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Custom Digital Filters for Dynamic Signal 
Analysis 

by Charles R, Panek and Steven F. Kator 



DIGITAL FILTERING offers some powerful advan- 
tages when used with a fast Fourier transform (FFT) , 
because it can perform frequency domain analysis 
in a narrow band around some arbitrary center frequency. 
To understand the need for a digital filter in dynamic signal 
analysis, let's first review the basics of FFT analysis and 
sampled systems, 

Basic FFT Analysis 

The fast Fourier transform is a special case of the discrete 
Fourier transform, a mathematical algorithm that translates 
a series of samples of a signal, evenly spaced in time, into 
a similar series of samples of the frequency spectrum of 
the signal. A finite set of these evenly spaced time samples 
is known as a time record. A time record with 1024 samples 
spaced 100 microseconds apart will be transformed into a 
frequency series or spectrum of 1024 evenly spaced sam- 
ples between dc and the 10- kHz sampling frequency f s . 
The frequency domain information for frequencies above 
iJ2 is redundant, that is, the spectrum is essentially folded 
around f s /2. This phenomenon determines the minimum 
rate at which a signal of a certain bandwidth can be sam- 
pled. Thus, provisions must be made to ensure that the 
sampled signal contains no frequency components of sig- 
nificance higher than $J2 so that the FFT results do not 
fold back on themselves, an effect called aliasing. 

A very basic FFT -based signal analyzer (Fig. 1) consists 
of a sampler preceded by a low-pass antialiasing filter with 
a cutoff frequency less than f s /2 f and an FFT processor that 
operates on fixed-length time records. But what if the user 
has an interest in signals at very low frequencies, or in 
signals spaced very close together in frequency? With a 
sampling rale as described above H this analyzer would only 
have about 10 Hz of frequency resolution. How can this be 
improved? One possibility is to lengthen the time record. 
A 10,000-point time record would allow T 1-Hz resolution. 
However, it would take 10 times as long to capture (he 
data, and roughly 13 times as long to calculate the FFT* 

Another possibility is to reduce the sampling frequency. 
This reduces the spacing of frequency components by a 
similar amount; a 1-kHz sampling frequency gives about 
1-Hz resolution with a 1024-point transform. This would 



be an ideal solution, except that the cutoff frequency of 
the low-pass input filter must then be changed to ensure 
that the spectral content of the sampled signal is entirely 
below half the new sampling frequency. To allow for many 
different band widths, we would have to be able to switch 
between many analog low-pass filters, or perform a compli- 
cated retuning of a single filter. 

There is another way to lower the sampling rate, how- 
ever, and that is by reducing the bandwidth of the signal 
after it is sampled and then discarding samples at regular 
intervals. This process is known as sample rate decimation, 
and the low-pass filtering is performed by a digital decima- 
tion filter, To increase the resolution by a factor of five, 
we would pass the output of our sampler through the digital 
filter, filtering the data to one-fifth the original bandwidth, 
and then discard four out of every five samples. We can 
do this because the filtered signal is now oversampled, that 
is, it contains redundant information. 

Digital filtering offers some other special advantages. The 
performance of a digital filter is a function of the coefficients 
used and the length of the data word. Once a digital filter 
design is completed, the bandwidth, flatness, noise, phase 
response, and other parameters of the filter are exactly 
repeatable from circuit to circuit. The adjustments and toler- 
ance variances characteristic of analog filters are eliminated. 

Zoom Operation 

By cascading decimating digital filters, we can reduce 
the bandwidth and increase the frequency resolution by 
larger amounts. We still have an important restriction, how- 
ever: the maximum frequency of our analysis continually 
drops lower. Suppose we wish to look at the modulation 
sidebands of a 4-kHz carrier frequency. If the modulation 
is very low in frequency, say 1 Hz, we will be unable to 
resolve the sidebands with the simple system of Fig. 1, 
Simply decimating the converter output by a factor of two 
will filter out the carrier, sidebands and all, assuming a 
10- kHz sampling frequency. 

The nature of the FFT limits us to frequency analysis of 
a time series whose frequency content extends from dc to 
half the sampling frequency. We can overcome this limita- 
tion by borrowing a technique from the designers of radio 
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Fig. 1, Block diagram of a basic 
fast Fourier transform signal ana- 
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receivers and swept -frequency spectrum analyzers. That 
is, demodulation or analysis is performed in a fixed inter- 
mediate frequency (IF) band, and any frequency band of 
interest is mixed with a local oscillator (LO) signal to trans- 
late it to the IF band. 

Here, we can consider our IF band to extend from dc up 
to half the sampling frequency, adjustable in bandwidth 
by the decimation rate of the digital filter. Thus, we can 
mix the 4-kRz carrier frequency with a 4-kHz local oscil- 
lator frequency and analyze the dc and low-frequency prod- 
ucts using a decimated sample rate that will give us the 
desired resolution, In addition, we can perform quadrature 
mixing, that is, multiply the input by sine and cosine waves 
at the LO frequency. By processing the resulting complex 
numbers in a complex FFT. frequency data both above and 
below the LO frequency can be obtained without generating 
unwanted images. We can once again use the advantages 
of digital signal processing to do this mixing in a stable 
and accurate way. The output of the analog-to-digital con- 
verter (ADC) is fed into a pair of digital multipliers, along 
with sine and cosine wave sequences of the desired center 
frequency sampled at the same rate. The resulting digital 
products are fed into a pair of identical digital filters for 
sample rate reduction, This is referred to as "zoom" oper- 
ation, for it allows the user to expand the view of a given 
portion of the input signal spectrum. 

Architecture and Algorithm 

We have arrived at the basic digital signal processing 
block diagram of the HP 3561 A Dynamic Signal Analyzer 
(Fig. 2). A 100-kHz analog antialiasing filter precedes a 
sample-and-hold circuit and ADC running at a sample rate 
of 256 kHz, This sample rate gives nicely spaced frequency 
points when used with 2 N -steed time records (in this case. 
1024 samples). The digital LO generates cosine and sine 
sequences, sampled at the ADQ sample rate, which are fed 
into the mixer multipliers of the real and imaginary com- 
ponent digital litters, respectively, along with the ADC out- 
put, In baseband mode t where no frequency translation is 



desired, the cosine output of the LO is a constant, and only 
the real component filter is used. The 20-bit output of the 
filter chips is rounded to 16 bits by the control 1C and fed 
into RAM, where it is accessed by the FFT processor. 

It was determined on the basis of noise performance and 
data scaling that the digital filters should use 20-bit num- 
bers to represent internal data. The A-to-D sample rate is 
specified at 256,000 samples per second. Hence, to allow 
serial processing of a decimated sample between each new 
ADC sample, a clock rate of 2 (samples) x 20 (bits) x 256 
kHz, or 10.24 MHz, would be required. 

Because of power limitations and other concerns, an al- 
ternative to strictly serial processing is used. Bits of the 
data word are processed in pairs of adjacent bits. Thus, the 
two leas! -significant bits jLSBs) are processed on the first 
clock cycle, the next two bits are processed on the second 
cycle, and so forth. This lowers the required clock rate to 
5.12 MHz. 

Mixer Multiplier 

The multiplier used in the HP 3561As signal processing 
circuitry uses a serial paired-bit shift-and-add scheme that 
operates on two ^-complement words. To understand its 
operation, consider a basic serial/parallel unsigned multi- 
plier as shown in Fig. 3. 

The circuit is initialized by clearing all the D flip-flops 
except the one flip-flop at the rounding point, which is 
set. and then loading the multiplicand into the parallel 
latch. As each multiplier data bit is brought in (least-signif- 
icant bit first], a new partial product is formed at the AND- 
gate outputs. The partial products are accumulated by the 
carry-save adders, and the result is presented serially at 
the output, with one new product bit for every new multi- 
plier bit, The circuit is clocked utili! all the product bits are 
shifted out, and the input data is padded with zeros as 
necessary, 

We may enhance this design to accommodate two Vcom- 
plement numbers by adding a serial negator ahead of the 
most-significant AND gate in the serial input path. If the 




256 kHz 



Fig. 2. Basic digital signal processing block diagram of the HP 356 1 A Dynamic Signal Analyzer 
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multiplicand is negative, this has the same effect as sub- 
tracting the last partial product. 

This design was adapted to work with bit pairs by split- 
the serial input data stream into odd and even bits, 
and generating three two-bit-wide data streams represent- 
ing the da la multiplied by 1 T 2, and 3. The parallel data 
(multiplicand) is also grouped into pairs. The AND gates of 
the basic circuit are replaced with pairs of 4:1 multiplexers, 
with the multiplicand bit-pairs driving the select lines and 
the x i f x 2, and X ij multiplier bit-pairs driving the inputs. 

The final design is shown in Fig, 4. The ADC and LO 
data words are brought in serially and converted to bit 
pairs. The LO word enters a two-bit- wide, serial -in, paral- 
lel-out shift register. The ADC word passes through a hard 
limiter that prevents overflowed ADC data from entering 
the filter, and a 0.75 scaler, which reduces the amplitude 
of the data so that overshoot in the filter will not result in 
overflows. This multiplier outputs the 20 most-significant 
[.Ms iMSBs) of the product of a 16-bit LO word and a 13-bit 
ADC data word in ten clock cycles, which suits the require- 
ments of the digital filter perfectly. 

Cascaded, Pipelined, Second-Order Sections 

There are two filters on the HP 3561 A" s custom digital 
filter chip. One, with a bandwidth of approximately 20% 
of the input sample rate, is used before reducing f 4 by a 
factor of 2. The other filler, with a 7S% bandwidth, is used 
before decimation by a factor of 5. Both of these filters are 
implemented as a cascade of second-order sections. Both 
filters also share the same digital circuitry. Filter selection 
is accomplished by an additional address bit to the on-chip 
ROMs that contain the filter coefficients, Each filter is an 
eighth-order elliptical design, with coefficients chosen to 
meet iJgsigjQ criteria while minimizing the number of hils 
to represent them. Having fewer bits in the coefficients 
reduces the amount of digital hardware needed to imple- 
ment the filters. 

Pipelining of the circuitry is used extensively between 
the second-order sections of the filter. That is r data is de- 
layed by one sample between inir.h filter section. Thus, the 
entire filter processes several samples at one time, which 
increases throughput without increasing clock frequency. 

There are ten channels of data storage associated with 
each filter section and the output of the filter, Using this 
memory, data can be recycled through the filter up in ten 
times after the sample rate has been reduced each time. By 
processing decimated data between samples from the ADC t 
the filter sample rate need only be twice that of the A-tu D 
converter 1 

Peled-Liu Second-Order Sections 

Each second-order filter section must implement a sim- 
ple difference equation: 



y n = a x n + a^ t + a 2 x n 



biYn-y " b a y ft _ 



(1) 



where x n is the nth input word, y n is the nth output word T 
and a^ and b, are the filter coefficients. 

To implement this equation with a minimum of 
hardware, a technique variously known as the Peled-Litr 
or Princeton multiplier is used- A brief explanation of this 



technique and its adaptation to this design follows. If we 
assume that data is represented by twos-complement bi- 
nary numbers of B bits bounded by 1 and -1. we can 
represent x a of Equation 1 by: 



where x/. is the jth bit of the nth input word, A similar 
expression represents y„ of Equation 1, Using these expres- 
sions, the difference equation can be rewritten as: 

|fxU2 '-xJU -b, ^yA-^-i-y^ 



+ a 



[3] 



We can define a function & with five binary arguments as: 
<//fx\ x 2 , x 3 . x 4 , X s ) = a tJ x 1 + a 1 x 2 + a 2 x 3 -b i x 4 -b 2 x 5 (4) 

where x k = or 1 
Then Equation 2 can be rewritten as: 

ii i 
5 = 1 

-.Hx^x^.x^.y^.ylU) (5) 

The function tjtcen have only 32 possible values, depend- 



Accu mutator 
and Shifter 




>YJU 



Fig, 5. Basic Peled-Ltu second -order filter. Blocks A. B. and 
Dare seriai-in, serial-out shift registers, Block Cis a parallel-in, 
serial -out shift register. 
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ing on its five binary arguments. Since the values of a t and 
bj are fixed in each filter section, this suggests implement- 
ing this function in ROM. The entire second-order section 
is shown in Fig, 5, 

Shift registers A through U are all B bits long. The nth 
input word x„ isfed T least-sign if ic:anl-bit first, into the ROM 
and shift register A, The ROM is simultaneously addressed 
by ihe corresponding bits of the x ri _ P x ri _^ y ri _j and y n _ 2 
words. The ip value generated on each clock cycle is ac- 
cumulated with the result of the previous operation, shifted 
right one hit. On the first cycle, the accumulator register 
is preloaded with a rounding value, which is added to 
assure minimum bias when the result is truncated back to 
B bits. On the final cycle of each filter operation, the if/ 
value corresponding to the sign bits of the data word is 
subtracted from the accumulator, and the accumulator 
value is parallel-loaded into shift register C 



Paired-Bit Implementation 

To maintain a lower clock rate, as mentioned earlier, the 
design of Fig. 5 was adapted to process hit pairs, using tw r o 
Princeton multipliers per second-order section. Equation 
5 can be rewritten: 

B 1 

i i.ri.... 

H- 2 

| J 4.. 



*#&4-l^ frj!U:^r*l 



(6) 



where R is an even number. 
This leads to a second-order section design as shown in 
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Fig. 6, Basic circuit of Fig. 5 
adapted to process bit pairs using 
two 10-channet memory banks 
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Fig. 6. Memory is organized as banks of 10-bit shift registers. 
One bank stores the even bits of the 20-bit data words, and 
the other stores the odd bits. Each bank consists of 10 
channels of memory, one for each pass the data can make 
Ihrough the filter. Thus, a second-order section can filter 
one sample from the outpnt of the mixer multiplier, and 
then store away its state and operate on a sample from a 
data sequence thai has already passed through the filter 
and had its sample rate reduced, 

Each Princeton multiplier generates a partial product, 
20 bits long, after rounding. These numbers are then added, 
two bits at a time, in the carry /save adder. The resulting 
sum is then spht into its odd and even bits, and fed into 
the appropriate channel of y n -^ memory. 

The y n _ T memory for the first filter section serves double 
duty as the x n _j memory for the second section. Pipeline 
registers are used for delaying and rcsynchronizing the data 
between sections. 

The simplified block diagram for the complete filter is 
given in Fig. 7 + Note that actually five second-order sections 
are used to implement an eighth-order filter. The first sec- 
tion implements a transfer function consisting of a scaling 
constant and two poles, and the last section implements a 
transfer function with only a pair of zeros. This implemen- 
tation was chosen to minimize circuit area on the !C. 

The last section of the filter is followed by an overflow 
detect and hard-limit circuit. Although the data is scaled 
down at the input so that internal overflows cannot occur, 
to preserve unity gain, the data is scaled up in the last 
section. If an overflow does occur there, the output data is 
set to the appropriate ± full-scale value, 

Once data emerges from the hard limiter, it is stored in 
the lD-channel-by-2 0-bit feedback memory » where it is re- 
circulated until new data is written into the same channels. 
Thus, sample-rate reduction of filtered data is accom- 
plished by merely writing into feedback memory more often 
than reading from H. 

Limit Cycles 

Just as there are challenges that face the designer of an 
analog filter who wishes to maximize noise performance 
and dynamic range, the digital filter designer must cope 
with problems unique lo digital signal processing. These 
include noise-like signals and dc offsets generated within 
the filter because of errors In rounding, and low-level self- 
oscillations in the absence of a signal. 

The phenomenon of limit cycle oscillations occurs in 
many digital signal processing systems where there is a 
rounding operation within a feedback loop. As a simple 
example, consider the following difference equation: 



x n + 0.9y n _ 1 



I7| 



If we set y n _ , = and x„ = 1 when n = 0. and x n = when 
n>0. then the following sequence results: 

Vo = 1 
y, = 0.9 
y 2 = 0.81 
y 3 = 0.729 
y 4 = 0.6561 



As can be seen, the sequence decays exponentially to- 
ward zero. Now let's limit the precision of the arithmetic, 
and round each result to the nearest tenth. Then: 



y 


= 1 


7t 


= 0,9-* 0.9 


Yz 


= 0.81 -* 0.8 


y* 


= 0.72 -*0.7 


y 4 


= 0.63 — 0.6 


y& 


= 0.54 — 0.5 


y B 


= 0.45 — 0.5 


¥? 


= 0.45 — 0.5 



This time the sequence "hangs up" at 0.5; even with no 
input I here is an output, A dc output like this is known as 
a deadband. When there are two feedback terms in our 
difference equation as in Equation 1 [thus T implementing 
a pole pair), the output can decay into an ac oscillation 
known as a zero-input limit cycle. Bounds for the 
amplitudes of these limit cycles have been developed. 3 In 
general, they relate to the Q of the poles of the transfer 
function. As the Q of the poles increases, limit cycles of 
higher amplitude become more likely, Since sharp-cutoff 
filters tend to have high-Q poles, a digital filter designer 
can expect to have problems with limit cycles when trying 
to maximize performance. 

When incorporating a digital filter in an FFT analyzer, 
the presence of limit-cycle oscillations is especially unde- 
sirable. Because these oscillations tend to be roughly 
sinusoidal near the resonant frequency of the pole pair, 
they will appear as a spurious signal in the frequency dis- 
play of the analyzer when the input is removed or held at 
very low levels, To add to the confusion, the oscillations 
can have different values or go away altogether, depending 
on the trajectory the filter input takes as it decays. 

Limit cycles are not readily predictable in complicated 
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Fig. 7, Simplified block diagram of 
complete digital filter, 
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fillers, For a basic second-order filler, it is a simple matter 
lo find all the zero-input limit cycles by repeatedly allowing 
the filter output to decay (or settle into oscillation) after 
presetting the filter's registers to every possible combina- 
tion of small values (one pair at a time). However, the 
results of such an experiment do not guarantee the limit 
cycle performance when the filler input is other than zero, 
or when the filter is cascaded with other second-order sec- 
tions, 

Deadbands and limit-cycle oscillations t^m be explained 
qualitatively as the rounding operation dominating the fil- 
ler coefficients in determining the decay of the filter's out- 
put. A design tactic that has often proven successful is to 
inject a low-level dither signal at the rounding point. In 
the example shown above, subtracting a value of 0.U5 from 
the number before rounding will allow the filter oulput to 
drop below 0,5. If we add 0.05 as often as we subtract it, 
the ro u tiding operation will still yield a minimum dc offset 

This solution has its disadvantages, however, [fa simple 
periodic dither signal is used, it too will appear as a spuri- 
ous signal on the output of the FFT stage, albeit of predict- 
able amplitude and frequency- As a matter of practice, such 
a dither signal often must be of an amplitude equal to or 
greater than the ampl itude of the oscillations it must break up. 

In the HP 3561 A Dynamic Signal Analyzer, limit cycles 
are eliminated by using a dither signal derived from a long 
pseudorandom noise generator, This proves quite effective 
in eliminating deadbands and limit cycles altogether. The 
energy from this noise generator is not concentrated at one 
frequency, but is spread evenly over the entire spectrum 
when its sequence length is sufficiently long. Thus, instead 
of a spurious component at some frequency on the FFT 
display, the dither appears as an increase in the noise floor. 
In fact, the peak value of the pseudorandom noise dither 
at each rounding point in the filter is less than one LSB of 
the 20-bit data word. Therefore, not only are there no limit 
cycles or deadbands, but the dither generator is not visible 
at the output under zero input conditions. 

Offsets 

Offsets or errors as a result of rounding can result in 
significant measurement inaccuracies. Offsets in the filter 
output are most problematic during zoom measurements, 
because input signals at the same frequency as the local 
oscillator are downconverted to dc and are displayed at 
midscreen. A dc error in the filter output thus results in a 
spurious component at the center frequency on the d is p 1 , i \ 

One would expect that since the filter data is 20 bits 
long, any reasonable rounding scheme should result in 
acceptably low distortion. However, several effects exacer- 
bate the problem. 

Since the HP 3561 A filter is a cascade of five sections, 
and since the decimation algorithm entails further cascad- 
ing the filter with itself up to ten times, any consistent 
offset in the section arithmetic will accumulate. Also, since 
rounding occurs after each multiply within the second- 
order sections, and because of the quantization of the coef- 
ficients, the results of certain coefficient multiplications 
have nonuniform statistical distributions. Hence, biases 
occur. Finally, the dc gain from some internal nodes to the 
output of the filter can be fairly high, particularly with the 



higher-Q, divide-by-5 coefficients, 

Curing DC Offsets 

As a first example, consider the problem of rounding 20 
bits to 16 bits. If we simply truncate the unwanted bits, 
the magnitude of the error introduced will be distributed 
between and - 15 LSBs of the 20-bil word. If we assume 
that the value of the truncated bits is uniformly distributed, 
the error will have a mean of - 7.5 LSBs (see Fig. 8a). 

The traditional solution is to add a value of 8 LSBs (VS» 
LSB of the 16-bit result), which yields a mean error of Vs. 
LSB of the 20-bit word. This is quite good, but not good 
enough for our application. 

Also consider the case of nonuniform! y distributed val- 
ues that are truncated. For example, suppose the first bit 
below the truncation point has a 75% probability of being 
zero. The truncation error will be distributed as in Fig. 8b, 
with a mean of - 5.5 LSBs. If we do the standard procedure 
of adding ft LSBs, an error of 2,5 LSBs still exists. 

The solution employed in the HP 35 61 A digital filter 
consists of adding one of I wo predetermined values before 
truncation, The two values are calculated so that their mean 
added to the mean truncation error equals zero* The choice 
;il which value to add is controlled by a pseudorandom 
register bit with a 50% duty cycle. Thus* on the average, 
the rounding error will be zero. This lechnique has two 
major advantages. It is flexible enough to deal with highly 
nonuniform distributions t and it allows the designer to 
adjust the variance in the rounding (noise power) to suit 
such needs as breaking up limit cycles. 

To see how this so-called "dithered rounding" is applied, 
return to the first example above (Fig. 8a). If we alternately 
add zero and 15 LSBs to the value before truncation, there 
will be an average error of zero in the operation. Likewise 
in the second example [Fig. 8b), we may alternately add 
zero and 11 LSBs before truncating, again achieving zero 
mean error. Clearly, this has profound effects on the round- 
ing noise, but the average noise power will generally be 
less than four times that of simple truncation, The use of 
random alternation is important, because the alternative 
use of periodically determined rounding is equivalent to 
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Fig. 8. Example of probability mass function of truncation error 

for (a) uniformly distributed values and(b) nonuniform^ distrib- 
uted values of LSBs when truncating a 20-bit word to 16 bits. 
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modulating the signal, which is not normally desirable in 
a digital filter. 

Mask Design Process 

The digital decimation filters used in the HP 3561A are 
implemented as a set of three custom NMOS integrated 
circuits, fabricated at HP's Luveland Technology Center. 
Two identical filter chips, containing 30,000 transistors 
each, perform simultaneous filtering of the real and imag- 
inary data. A control chip containing about 10.000 transis- 
tors provides correct sequencing of control signals for the 
filter ICs and provides an interface between them and the 
main instrument processor and the memory used by the 
FFT processor. 

Designing integrated circuits of this complexity requires 
a methodology that is logical and well thought out, and 
that Leaves no room for errors. Once the filter algorithms 
and topology are determined, and the digital design com- 
pleted, the process of translating schematics into silicon 
begins, This process starts with preliminary layouts of crit- 
ical portions of the circuitry, both to help prove feasibility 
of the design, and to gain insight into the area requirements 
of the circuit. A large portion of the circuitry in the digital 
filter consists of shift register memory, A simple shift regis- 
ter cell layout was optimized for size, and, along with a 
few other cells, this critical cell was used to develop an 
architectural plan showing the sizes and locations of the 
major blocks of circuitry and the power and signal routing 
between these blocks. 

The artwork for the chips is built in a hierarchical fash- 
ion. The chip Is divided into large functional blocks, (mul- 
tiplier, filter arithmetic, memory, et cetera). These blocks 
are rectangular and their horizontal and vertical dimen- 
sions match. The large blocks are themselves composed of 
smaller, regularly shaped blocks of medium complexity 
(control circuitry, parallel arithmetic, input latches, et cet- 
era), This hierarchy continues down to the most basic 
blocks of circuitry, called leaf cells (e.g.. shift registers and 



flip-flops). 

A keystone of this design method is the use of a standard 
metal pattern. The metal layer on the IC is used to distribute 
power and route signals over long distances. Power distri- 
bution on the chip is done using a metal pattern than looks 
like a pair of interdigitated combs (Fig. 9a), One set of 
fingers is connected to the 5 V supply, the other to ground. 
Thus, the ground and 5V lines are spaced at regular inter- 
vals, alternating across the chip. Another element of the 
standard metal pattern is that room for a few metal signal 
lines is left between the power lines. 

All leaf cells are built using the standard metal pattern 
(Fig. 9b). Every cell must be the same height and includes 
half of a 5V line at the top and bottom of the cell, and a 
ground line across the middle of the cell. Between I he 
power lines is room for up to four metal signal bus lines. 
Circuitry is laid out underneath these metal lines, using 
either the metal bus lines or portions of the underlying 
polysilicon Layer for interconnect. A half-height cell can 
also be used for smaller circuit blocks. Thus, chip area that 
might otherwise be used just for busing signals from one 
block to another can also contain circuitry. Care must be 
exercised in layout, of course, to ensure that noise-sensitive 
nodes are not placed under clock or data lines* 

Production Considerations 

It is not sufficient that an IC design be functional in an 
instrument under nominal conditions. Design goals for the 
chip in List be developed early that will permit reliable op- 
eration given variances in ambient tempera I ore, supply 
voltages, fabrication processing, and the performance of 
external circuitry, 

Testing. Given these criteria, production testing of these 
integrated circuits is a two-part process. Each IC must he 
thoroughly checked for correct functional operation of all 
of its many thousands of transistors, then this operation 
must be confirmed under test conditions that attempt to sim- 
ulate worst-case conditions of operation in the instrument. 







(a) 



fb) 



Fig, 9. (a) Simplified layout of top 

metal pattern for digital fitter chip. 
(b) Standard metal pattern used in 
designing cells for the digital filter 
chip (see shaded area in (a)} 
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Complete functional testing of large, complex chips pre- 
sents a special challenge, especially with a circuit like the 
digital filter, whose data inputs are all serial. Determining 
ei sequence of ones and /ems on these inputs that is guaran- 
teed to excite all of the many thousands of nodes oo the 
chip is a virtually impossible task, given the length of time 
a complete simulation of the entire chip would take. 

Because most nf the internal data path of the digital f liter 
is also serial, direct accessibility of internal nodes is not 
the problem it might otherwise be. If any node in the serial 
data path is stuck at a one or zero, this condition will 
propagate an error to the output. Other, nonserial portions 
of the circuit such as ROMs, parallel adders h and control 
circuits, are tested by operating the chip in each of its 
modes, and exciting the filter with signals of large dynamic 
range for as long a test as Is feasible. In addition, a special 
test mode allows use of any of the 10 channels of data 
memory without passing data through the previous chan- 
nels Efrst, thus speeding up the memory test and helping 
to isolate any problems to a specific channel. 
Test Program Generation. Rather than manually generate 
a set of stimulus and response vectors to test the chips and 
then code these vectors into test system programming Ian- 
guage one clock cycle at a time, a software system was 
developed that automatically generates test system Ian- 
guage code when given a stimulus sequence. This system 
generates response vectors using the same software 
simulator that was used to model the chips during the 
design process. Other software then converts these vectors 
into serial data streams used by the chips and generates 
digital test code for the wafer and packaged parts test sys- 
tem used in production. 



pool, reliable operation even under conditions of elevated 
ambient temperature. The heat sink mounts to the board 
using a spring loading scheme that keeps the heat sink's 
mating surface flush and in firm contact with the surface of 
the package without unduly stressing the brittle ceramic 
package. 

Acknowledgments 

Andy Purceli did the circuit design and much of the 
layout design of the control chip- Wayne Gravelle brought 
his many years of experience in layout to bear in executing 
the filter chip. His concepts of hierarchical design made 
an overwhelming task manageable. Dick Roberts of the Uni- 
versity of Colorado gave valuable assistance in the design 
of the filter architecture, Lynn Schmidt provided much 
guidance during the early stages of the project. Charlie 
Potter did the groundwork for the control chip and man- 
aged the development of both ICs, keeping the road ahead 
clear. 

Special thanks go to all the folks at the Loveland Tech- 
nology Center. In particular, thanks to Larry Chesler and 
Doug Brirtlett for their support and assistance during the 
design phase, to Wayne Gisel and Hal Cook for helping to 
bring these two designs into production, and last, but cer- 
tainly not least , kudos to Sandy Sherlzer and her crew in 
the computer-aided artwork department for their many 
long hours and patience with the last-minute changes of 
engineers. 

The HP Design Aids group, along w r ith the customer ser- 
vice group at the Loveland Technology Center, developed 
a set of software tools that helped to close the design loop 
on these complicated ICs + 



Packaging, Cooling 

Dense, high-speed chips like the digital filter and control 
chips require special consideration for packaging and cool- 
ing. A 64-pin T grid-array package was chosen for both parts. 
This choice not only provides the needed number of pins 
for the control chip* but offers a low r thermal resistance 
between the substrate and the top surface of the package. 

A cast aluminum heat sink that fits over the chips on the 
printed circuit board in the HP 3561 A Analyzer promotes 
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